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HALO GAS AND GALAXY DISK KINEMATICS DERIVED FROM OBSERVATIONS AND ACDM SIMULATIONS OF 
Mgll ABSORPTION SELECTED GALAXIES AT INTERMEDIATE REDSHIFT 
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ABSTRACT 

We obtained ESI/Keck rotation curves of 10 Mgll absorption selected galaxies (0.3 < z < 1.0) for which 
we have WFPC-2/HST images and high resolution HIRES/Keck and UVES/VLT quasar spectra of the Mg II 
absorption profiles. We perform a kinematic comparison of these galaxies and their associated halo Mgll 
absorption. For all 10 galaxies, the majority of the absorption velocities lie in the range of the observed galaxy 
rotation velocities. In 7/10 cases, the absorption velocities reside fully to one side of the galaxy systemic 
velocity and usually align with one arm of the rotation curve. In all cases, a constant rotating thick-disk 
model poorly reproduces the full spread of observed Mgll absorption velocities when reasonably realistic 
parameters are employed. In 2/10 cases, the galaxy kinematics, star formation surface densities, and absorption 
kinematics have a resemblance to those of high redshift galaxies showing strong outflows. We find that Mgll 
absorption velocity spread and optical depth distribution may be dependent on galaxy inclination. To further 
aid in the spatial-kinematic relationships of the data, we apply quasar absorption line techniques to a galaxy 
(v'c =180 km/s) embedded in ACDM simulations. In the simulations, Mgll absorption selects metal enriched 
"halo" gas out to ^ 100 kpc from the galaxy, tidal streams, filaments, and small satellite galaxies. Within the 
limitations inherent in the simulations, the majority of the simulated Mg II absorption arises in the filaments and 
tidal streams and is infalling towards the galaxy with velocities between -200 < Vr < -180 km s"'. The Mgll 
absorption velocity offset distribution (relative to the simulated galaxy) spans 200 km s"' with the lowest 
frequency of detecting Mg II at the galaxy systematic velocity. 

Subject headings: galaxies: halos — galaxies: kinematics and dynamics — galaxies: intergalactic medium — 
quasars: absorption lines 



1. INTRODUCTION 

In a cosmological context, galaxy formation occurs via ac- 
cretion of gas from the cosmic web and from galaxy-galaxy 
mergers. The galaxy kinematics reflect these processes. As 
galaxies evolve, a complex interplay develops between the 
star formation, which can generate winds, and an array of 
kinematic structures such as tidal streams, galactic fountains, 
and filamentary infall that comprise an extended gaseous halo. 
This overall picture is suggested by observations and ACDM 
cosmological simulations of galaxy formation. However, we 
lack a thorough understanding, both observationally and theo- 
retically, of how these processes precisely affect the dynamics 
of galaxies and their extended halos. High quality, detailed 
observations are required to further develop this working sce- 
nario and produce a comprehensive model of galaxy evolution 
in the cosmological context. 

Observations of local galaxies provide detai led views of gas 
disks a nd the inner ~ 25 kpc of their halos. lOosterloo et al.l 
(|2007) obtained deep H I observations of NGC 891 and found 
lagging, differentially rotating halo gas kinematics with ve- 
locities decreasing with distance above the galaxy plane. 
Lagging halos are observ ed in several other local galaxies 
(e.g" LS ancisi et al. 2001; ISwaters et all Il997t iRandl 120001: 
iHeald et al.i i2007i) . There are also cases where halo gas 
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has been detected with velocities opposite to galaxy rotation 
(Fraternali et al. 2001; Oosterloo et al. 2007). 

Evidence for the accretion of cold gas, which may lower 
angular momentum and play a role in the development of 
lagging halos, has also been rapidly accumulating. Several 
galaxies are observed to have galactic fountains, and be sur- 
rounded by H I clou d complexes, mino r merger tida l tails, and 
IGM filaments (e.g. jHeald et"aLll2007l: ISancisi et al.ll2008h . 

Quasar absorption lines provide powerful probes of halo 
gas kinem atics to l arge galactocentric distances. In three 
galaxies. Cote et al.l (12005 ) found that low column density 
Lya absorption does not follow galaxy rotation to D = 
390 k pc. They suggested that the gas arises from the cosmic 



web. Bowen et al.l (l2002h found that Ly a absorption strength 
correlates with galaxy over density. At large distances, the 
halo gas kinematics may no longer be coupled to the galaxy 
kinematics, but reflect the motions of cosmic web. 
Quasar absorption line studies of the Mgll 

AA2796,2803 doublet produc ed by gaseous h alos of 

foreground gala xies (e.g., iBergeron & Boissl 119911: 
iLe B run et aP 1 1 993t LS teidel, Dick inson. & PerssonI 



1994; IChurchin. Steid el. & Vogt 

Churchill, Kacprzak. & Steidell 120051: fZibetti et alj 
Kacprzak et al. .2007t) probe the low ionization metal en- 
riched gas-galaxy dynamics and IGM interface. With 
Mgll absorption lines, we can study the kinematic con- 
ditions of galactic halos over a wide range of redshifts 
out to projected galactocentric radii of several hundred 
kpc. Mgll absorption arises in low ionization, metal 
enriched gas with neutral hydrogen column densities of 
10'^ < N(Hl) < 10^^ cm"^, and thus selects a large dynamic 
range of gas structures in the environments associated with 
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galax ies dRigbv. Charlton. & Churc"hnill2002t IChurchiU et al] 

mmh . 

The idea that Mgll absorbers could arise from spherical 
infall, disk-like rotation, or both, is also a topic of much 
interest. In a s mall sample of high reso lution Mgll ab- 
sorption profiles, iLanzetta & Bowenld 19921) inferred that ro- 
tation kinematics dominated at smaller impact parameters, 
whereas infall kinematics dominated with increasing impact 
parameter. Armed with a larger high resolution sample, 
ICharlton & ChurchiTll ( Il998h applied statistical tests to a vari- 
ety of kinematic models and concluded that pure disk rotation 
and pure halo infall models are ruled out. However, models 
with contributions from both disk rotation and spherical in- 
fall statistically reproduced absorption profiles consistent with 
observed kinematics. 

Among the more extreme structures, are galactic winds 
generated from star forming galaxies. At z < 1.4, 
iTremonti etsU ( l2007h detected Mgll blueshifted 500- 
2000 km s~' relati ve to post-starburst host galaxies and 
IWeiner et al.l (l2009l) found 300- 1000 km s"' blueshifts in 
Z ~ 1.4 star forming galaxies. In other ionic species, simi- 
lar outflows have been observed at z ^ 3 (Pettini et al. 2001; 
Shanlev et al.' '2003': LSteidel etal.ll200l iSimcoe et al.li2006i: 
Cabana c et al. 2008). None of these surveys have studied 
the dynamics of the galaxies themselves, which are clearly 
needed in order to obtain a full picture of the galaxy-halo dy- 
namics. 

A direct comparison of the galaxy disk kinematics and 
absorbing Mgll halo gas kinematics h as been performed 
for six z ^ 0.6 hig h ly inclined galaxies (ISteidel et al.l 120021: 
lElhson et al.1l2003l) . lElUson et al] (|2003|) found that the sys- 
temic velocity of a galaxy coincided with the center of the 
absorption system, which spanned more than 100 km s"' 
about the systemic velocity. Bond et al. (2001) used expand- 
ing shell models to explain that this peculiar absorption pro- 
file is lik ely caused by expan ding supernovae-driven super- 
bubbles. ISteidel et atl (l2002l) found that, in four of the five 
cases, the velocities of all of the absorption components lie to 
one side of the galaxy systemic redshift. The fifth case had 
a narrow, weak absorption centered at galaxy systemic veloc- 
ity. Since the halo gas velocities align in the same sense as the 
galaxy rotation, the velocity offsets of the absorbing gas rel- 
ative to the galaxy systemic velocity strongly suggest "disk- 
like" rotation of the halo gas. Using simple disk halo models, 
ISteidel et al.l (120021) concluded that an extension of the disk 
rotation with a lagging halo component (based upon proper- 
ties of local galaxies' halo gas kinematics) was able to explain 
some of the gas kinematics. However, the models were not 
able to account for the full velocity spreads of the gas. 

From a theoretical stand point, semi-analytical 

models and isolated ga laxy simulations (e.g.. 

Mo & Miralda-Escu dg 119961; ISur kert & Lin 2000; 
Lin & Murray 2000; Mailer & Bullock 2004; Ch en & Tinken 
i308; Kaufmann et al. 2008; Tinker & Chen 2008) have been 
invoked to study isolated galaxy halos. In these models, 
Mgll absorption arises from condensed, infalling, pressure 
confined gas clouds within the cooling radius of a hot halo. 
These models are quite successful at reproducing the general 
statistical properties of the absorber population. However, 
they lack the important dynamic influences of the cosmic 
structure and local environments. 

ACDM simulations have been able to synthesize the forma- 
tion a nd evolution galaxies within large scale structures. Re- 
cently, |Ceverino]&^i^2^(|2QQ9t), were able to naturally cre- 



ate, without ad-hoc recipes, extended galactic scale outflows 
and metal enriched multi-phased galactic gas halos. This was 
accomplished by studying the detailed physics of the forma- 
tion and evolution of the multi-phase ISM in parsec resolution 
simulations. These same prescriptions were then successfully 
applied in their large scale cosmological simulations. 

Since these cosmological simulations include all the poten- 
tial structures that can influence halo gas dynamics and in- 
clude the local environment, they provide a promising tech- 
nique for understanding the role of gas in galaxy evolution. 
The quasar absorption line method can be applied to simula- 
tions to examine structures selected by species such as Mg II, 
in the vicinity of galaxies. The goal is to compare directly 
observed absorbing halo gas kinematics and host galaxy kine- 
matics to those extracted from the simulations. In order to 
arrive at a deeper understanding, the observations should tar- 
get redshifts where detailed high quality kinematics can be 
obtained for a sample of galaxies with a wide range of orien- 
tations with respect to the quasar line of sight. 

We have obtained ESI/Keck rotation curves of 10 interme- 
diate redshift (0.3 < z < 1 .0) galaxies for which we have high 
resolution HIRES/Keck or UVES/VLT quasar absorption pro- 
files of Mgll, as well as WFPC-2/HST images. In this paper 
we perform a kinematic comparison of 10 galaxies and their 
associated halo Mgll absorption. We define halo gas to be 
metal enriched structures that give rise to Mgll absorption 
such as extraplanar gas, outflows, tidal streams, filaments, and 
satellite galaxies. We compare our observations with a simple 
rotating thick disk halo model (similar to the one employed 
by Steidel et al. 2002) and w ith the cosmological simulations 
of Ceverino & KlypinI (|2009|) . 

The paper is organized as follows: In § |2] we present our 
sample, and explain the data reduction and analysis. In §[3] we 
present the results of our galaxy-Mg II absorption kinematic 
observations, and in § |4l we compare the observed absorp- 
tion velocities with a simple disk kinematic halo model. In 
§ |5] we discuss the details of the cosmological simulations. 
We study a simulated galaxy and its halo structures in detail. 
We analyze the integrated total hydrogen and Hi column den- 
sity maps and the absorbing gas velocity distributions. We 
use these results to infer possible structures and kinematics 
drivers of that produce the observed Mg II absorption profiles. 
We also compute the star formation rate and star formation 
surface density of the simulated galaxy and compare them to 
previous observational results, and with two of the galaxies in 
our sample. We end with our conclusions in § |2l Throughout 
we adopt ah = 0.70, JIm = 0.3, flA = 0.7 cosmology. 

2. DATA AND ANALYSIS 

2.1. Sample Selection 

The selection of the sample presented in the study is 
based upon three steps. (1) We compiled a list of Mgll 
absorbers in high resolution (R = 45,000,FWHM ~ 6 km s"\ 
wifli £W,(2796) > 0. 02 A) H IRES/Keck (IVogt et alJ 
IT994h or UVESA^LT dPekker et al. 200Q) quasar spec- 
tra. We make no cut to the sample based upon equiv- 
alent width. (2) We then compiled all subsequent deep 
ground based imaging and spectroscopic redshift sur- 
veys of the quasar fields and selected the galaxies that 
have confirmed redshifts aligned with Mgll absorption 
(Bergeron & Bo issi 1 199 It Bergeron. Cristiani, & Shav^ 
Il992£ Steidel, Dickinson, & Persson 1994; Lowenthal et al] 
Il995t .Guillemin & Bergeron 1997; Chen et al. 19981 



Galaxy Disk and Halo Gas Kinematics: Observations and Simulations 3 



iLane etalJ [19981 This paper). (3) Finally, we selected 
galaxies for which WFPC-2/HST images were available 
and from which we ca n extract detailed g alaxy morpho- 
logical parameters (see Kacprzak et alj 12007) . The final 
sample of ten galaxies have an impact parameter range of 
26<D< 108 kpc. 

Our goal is to study the relationship between the spatial and 
kinematics relations between a galaxy and the Mg II absorbing 
gas in its vicinity. In our sample, two of the Mg II absorbers 
appear to be associated with galaxies of similar luminosity 
that exhibit signs of interaction in the HST images. These 
galaxies are Gl and G2 in the field of Q0450-132 and Gl and 
G2 in the field of Q1127-145 (see Figures |3fc,/7 and|5l re- 
spectively). The Q0450-132 galaxies show tidal asymmetries 
and have projected separated of 12 kpc. There are no addi- 
tional candidate galaxies in the image within 100 kpc (pro- 
jected) of the quasar. The Ql 127-145 galaxies also exhibit 
tidal asymmetries. Their projected separation is roughly 50 
kpc. In the Ql 127-145 field there is a third galaxy, G4, with 
much smaller luminosity, with a redshift that places it within 
70kms-i ofG2 and 25 km s"' ofGl. 

These particular systems, which are characterized by two 
roughly comparable luminosity galaxies with signs of inter- 
acting, pose an interesting challenge. They indicate that some 
Mg II absorption is arising in the complex environment of a 
major-major galaxy interaction. It is probably a fair statement 
to assert that interacting galaxies of roughly equal luminosity 
(and that is an important point) will have local environments 
very different than those of galaxies that clearly have no com- 
panion of comparable luminosity. It becomes an intractable 
problem to discern what portion of the absorption may be 
arising with gas associated with one or the other galaxy in 
such a pair. In cases where a single galaxy candidate can be 
assigned as the luminous host of the Mgll, it is possible to 
unambiguously study the spatial and kinematic relationships. 
Of course, it is always possible a very low luminosity coun- 
terpart is below the detection of the images; but such a com- 
panion would indicate a minor-major interaction and not a 
major-major interaction like the Q0450-132 and Ql 127-145 
pairs. Minor-major interactions would be more akin to the 
Magellanic galaxies in the 50 kpc vicinity of the Milky Way; 
they can be considered part of the Milky Way halo. Such a 
distinction would equally apply to the single galaxies in this 
sample if there is an unseen minor companion. 

In order to keep our sample as uniform as possible for 
the study of the spatial and kinematic connections between 
a galaxy and its associated absorbing gas, we limit our anal- 
ysis to the single galaxy candidates when ambiguities would 
arise, such as comparison with galaxy inclination, position 
angle, and impact parameter There is a very different nature 
to the major-major interacting pairs toward Q0450-132 and 
Ql 127-145 in that they may have a common gas envelope 
and therefore provide a slightly different probe of absorption 
galaxy properties. In cases where ambiguities do not arise in 
the analysis, we include all the galaxies in our sample. 

2.2. Quasar Spectroscopy 

Details of the HIRES/Keck and UVES/VLT quasar obser- 
vations are presented in Table [1] The HIRES spectra (except 
for Q0836 + 1 13) were reduced using IRAF''. The spectrum of 

' IRAF is written and supported by the IRAF programming group at 
the National Optical Astronomy Observatories (NOAO) in Tucson, Arizona. 
NOAO is operated by the Association of Universities for Research in Astron- 



TABLE 1 

Keck + VLT Quasar Observations 



Exposure 



QSO Field 


Zem 


Instrument 


Date (UT) 


(sec.) 


Q0002 + 051 


1.90 


HIRES 


1994 Jul. 05 


2700 


Q0229 + I3I 


2.06 


HIRES 


1999 Feb. 08 


3600 


Q0450-I32 


2.25 


HIRES 


1995 Jan. 24 


5400 


Q0454-220 


0.53 


HIRES 


1995 Jan. 22 


5400 


Q0836+1I3'' 


2.70 


HIRES 


1998 Feb. 26 


5400 


Q1127-145 


1.18 


UVES 


. . . <i 


24,900 


Q2206-199 


2.56 


UVES 


a 


53,503 



^ The Q 11 27- 145 quasar spectrum was obtained over multiple nights. The 
PIDs for this quasar are 67.A-()567(A) and 69.A-0371(A). The Q2206- 199 
quasar spectrum was also obtained over multiple nights for the following 
PIDs 65.0-OI58(A), 072.A-0346(A), and 074. A-020I(A). 
^ Data provided by Jason X. Prochaska jProchaska et aU2007^ ■ 



TABLE 2 
WFPC-2/HST OBSERVATIONS 







Exposure 




QSO Field 


Filter 


(sec.) 


PID/PI 


Q0002 + 051 


F702W 


4600 


5984/Steidel 


Q0229+131 


F702W 


5000 


6557/Steidel 


Q0450-132 


F702W 


2500 


5984/Steidel 


Q0454-220 


F702W 


1200 


5098/Burbidge 


Q0836+113 


F702W 


5000 


6557/Steidel 


Q1127-145 


F814W 


4400 


9173/Bechtold 


Q2206-199 


F702W 


5000 


6557/Steidel 



Q0836+ 113 was reduced using the MAKEE^ package. The 
UVES spectra were reduced using the standard ESO pipeline 
and a custom code called the UVES Po st-Pipeline Echelle 
Reduction (UVES POPLER, 'Murphy"2006). The quasar spec- 
tra are both vacuum and heliocentric velocity corrected. Anal- 
ysis of the Mgll absorption profiles was performed using 
graphic-based interactive s oftware of our own de sign (see 
IChurchiU et aWl99^ IIM IChurchiU & Yo^MMi for local 
continuum fitting, objective feature identification, and mea- 
suring absorption properties. The redshift for each Mgll 
system is computed from the opti cal depth weighted mean 
of the absorption profile (see .Churchill & Vog t ,2001) . The 
typical absorption redshift uncertainty is ~ 0.3 km s~'. The 
Mgll A2796 re st-frame equivalent widths are adopted from 
iKacprzak et al.l (l2C)07). Velocity widths of absorption systems 
are measured between the pixels where the equivalent width 
per resolution elemen t recovers to the 1 a detection threshold 
(IChurchiU etanfT999l) . 

2.3. HST Imaging and Galaxy Properties 

All WFPC-2/HST images were reduced using the WFPC- 
2 Associations Science Products Pipeline (WASPP^). The 
WFPC-2 astromet ry is calibrated to the USN02 Catalog 
(iMonet et al.lll998h . WASPP data quality verifications in- 
clude photometric and astrometric accuracy and correctly set 
zero-points. Details of the WFPC-2/HST observations are 
presented in Table |2] Galaxy apparent and absolut e magni- 
tudes are adopted from IKacprzak et al.l (l2007l l2008l) . respec- 
tively. The mfj()2w and mf-gMw magnitudes are based upon 

omy (AURA), Inc. under cooperative agreement with the National Science 
Foundation. 

* http://spider.ipac.caltech.edu/staff/tab/makee 

^ Developed by the Canadian Astronomy Data Centre (CADC) 
and the Space Telescope-European Coordinating Facility (ST-ECF): 
http://archive. stsci. edu/list/wfpc2/pipeline.html 
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TABLE 3 
Keck-II/ESI Observations 









Exposure 


Slit 


QSO Field 


^abs 


Date (UT) 


(sec.) 


PA 


Q0002 + 051 Gl 


0.851407 


2001 Oct. 16 


7200 


-8.5 


Q0229+131 Gl 


0.417337 


2006 Dec. 24 


6500 


134 


Q0450-132 Gl 


0.493937 


2006 Dec. 24 


5300 


-13 


Q0450-132 G2 


0.493937 


2006 Dec. 24 


5300 


-13 


Q0454 - 220G1 


0.483338 


2006 Dec. 24 


4800 


276 


Q0836+113 Gl 


0.786726 


2006 Dec. 24 


5300 


130 


Q0836+113 G2 




2006 Dec. 24 


5300 


130 


Ql 127 -145 Gl 


0.312710 


2006 Dec. 24 


3900 


129 


Ql 127 -145 G2 


0.312710 


2006 Dec. 24 


4200 


87.5 


Ql 127 -145 G3 


0.328266 


2006 Dec. 24 


600 


87.5 


Q2206-199G1 


1.017040 


2001 Oct. 16 


1800 


75 



There is no Mgll absorption associated witli Q0836+ 1 13 G2 in the litera- 
ture. Our HIRES data do not provide the necessary wavelength coverage. 



a) [Oil] Doublet 3730. 3727 For Q0454-220 zgal=0.4838 
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Fig. 1. — (a) A 2D spectral region around the [Oil] doublet from the 
z = 0.4838 galaxy in the quasar field Q0454-220. The [Oil] is spatially 
resolved and extends roughly 5" on the sky. The two sky lines extend the 
length of the 20" slit. — (b) A ID extraction of the above 2D image (thick 
dashed line in top panel) summed over 3 pixels in the spatial direction. The 
galaxy continuum fit is indicated by the dashed line. The 1 a uncertainty in 
the continuum is shown by the dotted lines bracketing the continuum. The 
sky lines and sky signal have been subtracted out. The solid line shows the 
Gaussian fit to the emission lines and the tick marks indicates the centroids 
of the line. 

the Vega syste m. As described in iKacprzalc et alj ( 12007b . we 
used GI1VI2D (ISimard et al.ll2002l) to model the galaxy mor- 
phologies, and measured the quasar-galaxy impact parame- 
ters, galaxy sky orientations, inclination angles (/), and posi- 
tion angles (PA) of their major axes with respect to the quasar 
line of sight. We fit each galaxy surface brightness profile 
with a Sersic bulge component (for 0.2 <n< 4.0) and an ex- 
ponential disk component. Additional modeled galaxy mor- 
phological parameters will be presented elsewhere (Kacprzak 
et al. 2010, in preparation). 

2.4. Galaxy Spectroscopy 



TABLE 4 

IVIgii Absorption And Galaxy Redshifts 



QSO Field 


Zahs 


Zgal 




Aw" 
(km/s) 


Q0002 + 051 Gl 


0.851407 


0.85180 


0.000066 


-66 


Q0229+131 Gl 


0.417337 


0.4167 


0.00020 


+135 


Q0450-132G1 


0.493937 


0.4941 


0.00015 


-33 


Q0450-132G2 


0.493937 


0.4931 


0.00012 


+168 


Q0454-220G1 


0.483338 


0.48382 


0.000066 


-98 


Q0836+113 Gl 


0.786726 


0.78682 


0.000028 


-16 


Q1127-145 Gl 


0.312710 


0.3132 


0.00020 


-112 


Ql 127 -145 G2 


0.312710 


0.3124 


0.00013 


+71 


Ql 127 -145 G3 


0.328266 


0.32847 


0.000027 


-46 


Q2206-199G1 


1.017040 


1.01655 


0.000013 


+73 



" A IV is the rest-frame velocity offset between the mean Mgll 
A2976 absorption line and the galaxy where, Aiv = c{zabs ~ 
Zga/)/(l+Zgfl/)kms"'- 



The ESI/Keck dSheinis et alj|2002h galaxy spectra were ob- 
tained over two nights; two were obtained in October 2001 
and eight were obtained in December 2006. Details of the 
ESI/Keck observations are presented in Table [3] For each 
galaxy, the slit position angle was chosen to lie along the 
galaxy major axis (except for Q450- 132 G2). The slit length 
is 20". Thus in some cases, we were able to simultaneously 
position two galaxies on a slit. Exposure times range between 
600 -7200s per galaxy. The wavelength coverage of ESI is 
4000 to 10,000 A, which allows us to obtain multiple emission 
Hnes (such as [Oil] doublet, H/3, [Olll] doublet. Ha, [Nil] 
doublet, etc.) with a velocity resolution of 1 1 km s"' pixel"' 
(FWHJVI - 45 km/s). 

In October 2001, the data were obtained with a 1" slit and 
2x2 binning. The mean seeing was 0.7" (FWHM) with par- 
tial cloud coverage. In December 2006, the data were ob- 
tained using a 0.75" slit with 2x1 binning. Binning by two in 
the spatial directions results in pixel sizes of 0.27-0.34" over 
the orders of interest. The mean seeing was 0.8" (FWHIVI) 
with clear skies. 

The data were reduced using the standard echelle package 
in IRAF. We used internal quartz illumination flat fields to 
eliminate pixel to pixel variations. In the science frames, sky 
subtraction was performed by fitting a polynomial function to 
each spatial column. A quasar or bright star spectrum in the 
same field was obtained and used as a trace in order to facil- 
itate the extraction of the galaxy spectrum. A spatially inte- 
grated spectrum was extracted in order to obtain an accurate 
galaxy redshift from the centroids of multiple emission lines. 
The lines listed in the legends of Figures|2}{6]were used to de- 
termine the galaxy redshifts. The rest-frame vacuum wave- 
length used for each emission line was obtained from the Na- 
tional Institute of Standards and Technology (NIST) database. 
Each spectrum was wavelength calibrated using CuArXe arc 
line lamps. Spectra were calibrated in IRAF using standard 
stars taken during the night of the observation. The flux is ac- 
curate to ^10% and we have made no corrections for slit loss, 
or Galactic reddening. 

In Figure [T] we show an example of a 2D spectrum of the 
spatially resolved [Oil] doublet from the z = 0.4838 galaxy in 
the Q0454 - 220 field. W e used a Gaussian fitting algorithm 
(see lChurchill et al.l l2000). which computes best fit Gaussian 
amplitudes, widths, and centers (redshift), to the galaxy emis- 
sion and absorption lines. The galaxy redshift was computed 
from the mean redshift of all the detected lines. Emission 
lines and absorption lines used to calculate the galaxy red- 
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TABLE 5 

Mgii Absorption And Galaxy Redshift Field Survey 





Galaxy 










W,(2T)f)) A 


Galaxies in 


QSO Field 


ID 


Zgal 


RctcrcncG 


D (kpc) 


Zahs 


This Study 


Q0002+05I 


Gl 


0.85180 


1 


25.9 ±0.5 


0.851407 


1.119±0.013 


X 




G2 


0.592 


2 


36.2 ±0.4 


0.591365 


0.102 ±0.002 






G3 


0.298 


2 


59.3 ±0.3 


0.298059 


0.246 ± 0.004 




Q0229+131 


Gl 


0.4167 


1,3 


37.5 ±0.5 


0.417337 


0.816 ±0.022 


X 


Q0450-132 


Gl 


0.4941 


1,2 


50.1 ±0.4 


0.493937 


0.674 ± 0.026 


X 




G2 


0.4931 


1,2 


62.7 ±0.7 


0.493937 


0.674 ±0.026 


X 


Q0454-220 


Gl 


0.48382 


1,4 


107.9 ±0.8 


0.483338 


0.426 ± 0.007 


X 


Q0454-220 


G2 


0.3818 


4 


103.4 ±0.3 




< 0.02 (3 cr) 




Q0836+I13 


Gl 


0.78682 


1,5 


26.9 ±0.9 


0.786726 


2. 148 ±0.023 


X 




G2 


0.48288 


1 


29.1 ±0.3 




. . . b 




Q1127-I45 


Gl 


0.3132 


1,3 


45.6 ±0.3 


0.312710 


1.773 ±0.006 


X 




G2 


0.3124 


1,3 


81.0±0.3 


0.312710 


1.773 ±0.006 


X 




G3 


0.32847 


1 


91.4±0.2 


0.328266 


0.029 ±0.003 


X 




G4 


0.3121 


6 


18.2±0.3 


0.312710 


1.773 ±0.006 




Q2206-199 


Gl 


1.01655 


1,7 


104.6 ± 1.4 


1.017040 


1.057 ±0.005 


X 




G2 


0.948 


7 


87.2 ±0.5 


0.948361 


0.253 ±0.002 






G3'^ 


0.755 


8 


44.2 ±0.7 


0.751923 


0.886 ±0.003 





" Galaxy Identification: (I) This paper , (2) 'Steidel, Dickinson, & Pensson' !l994}. (3) Berseron & Boisse j I991h . 
(4) Chenetal. a99 8j, (5) rLowenthal et aH fM s). (6) Lane et al. (1998), (7) Beraeron, Cristiani, & Shaver (199^, and 
f8) lGuiIIemin & BergeronI U997t ). We list the redshift for galaxies that were derived for this work. 

^ There is no Mgll absorption associated with Q0836+ 113 02 in the literature. Our HIRES data do not provide the necessary 
wavelength coverage. 

^ G3 was reported as a galaxy bv lGuillemin & BergeronI 119971) . However, our spectroscopic observations reveal that this 
um'csolved object is a Galactic star 



shift must have been detected at the 3 a level. The adopted 
redshift uncertainty for each galaxy was computed from the 
standard deviation in the redshifts computed from each emis- 
sion line. The galaxy redshifts are listed in Table IH their ac- 
curacy ranges from 2-45 km s~' . The galaxy velocity offsets 
from the optical depth weighted mean Mgll absorption are 
also listed in Tableland range from -12 to +168 km s"' . 

The rotation curve extraction was performed following the 
methods of Vogt et al. ( 1996) and Steidel et al. (2002). We ex- 
tract individual one-dimensional spectra by summing three- 
pixel wide apertures (corresponding to approximately one res- 
olution element of 0.82- 1.01") at one pixel spatial incre- 
ments along the slit. An error spectrum is also extracted for 
each of these apertures. To obtain accurate wavelength cali- 
brations, we extracted spectra of CuArXe arc line lamps at the 
same spatial pixels as the extracted galaxy spectra. Fitted arc 
lamp exposures (CuArXe) provided a dispersion solution ac- 
curate to ~ 0.15 A, or about 6.5 km s~' at the wavelengths of 
interest. Galaxy spectra are both vacuum and heliocentric ve- 
locity corrected for comparison with the absorption-line kine- 
matics. Each galaxy emission line (or absorption line in some 
cases) was fit with a single Gaussian (except the [Oil] dou- 
blet was fit with a double Gaussian) in order to extract the 
wavelength centroid for each emission line. 

An example of a three-pixel wide spectral extraction from 
the 2D spectrum is shown in Figure[TJ' where the [O ll] A3727 
and A3730 lines are detected at the 22cr and lla level, respec- 
tively (the significance level is the ratio of the measured equiv- 
alent width to the uncertainty in the equivalent width based 
upon error propagation using the error spectrum extracted for 
the same three-pixel aperture). The dashed line in Figure [T}? 
provides the spatial cut for which the spectrum in Figure ^ 
is illustrated. The velocity offsets for each emission line in 
each extraction were computed with respect to the redshift 
zero point determined for the galaxy (Table HJi. The rotation 
curves for the 10 galaxies obtained with ESI/Keck are pre- 
sented in Figures IIHS] 



3. DISCUSSION OF INDIVIDUAL FIELDS 

Here we discuss the halo gas and galaxy kinematics of ten 
galaxies in seven different quasar fields. In Table |5] we list 
all the galaxies in each field that have spectroscopically con- 
firmed redshifts. The table columns are (1) the quasar field, 
(2) the galaxy ID, (3) the galaxy redshift, (4) the reference(s) 
for the galaxy identification, (5) the quasar-galaxy impact 
parameter, D, and uncertainty, (6) the Mg II absorption red- 
shift, and (7) the rest-frame Mgll A2796 equivalent width, 
W, (2796), and uncertainty. Three new galaxies (Q0002 + 05 1 
Gl, Q0836+133 G2, and Ql 127- 145 G3) have been spec- 
troscopically identified in this work. 

Most of the quasar fields listed in Table |5] have been spec- 
troscopically surveyed for all galaxies with nir < 23, which 
translates to U ^ 0.08L,, fo r z - -5, out to ^ 1' (e.g. 
'Steidel. Dickinson. & Perss"onl Il994t iGuillemin & BergeronI 
1997; Chen et al. 199^. These works have been instrumen- 
tal in developing our current picture of the galaxy-absorber 
connection at intermediate redshifts. It remains possible that 
additional galaxies below this luminosity lurk in front of the 
quasar and may be associated with the absorbing gas. It is also 
possible that a less than 100% completeness in the confirma- 
tion of galaxy redshifts may result in an additional galaxy or 
galaxies also associated with the absorption. The concerns 
and caveats associated with incorrect idenitifications and the 
conclusion drawn for those works therefore also apply to this 
study. Galaxies fainter than nir ^ 23, or low mass galaxies 
hidden by the quasar PSF, such as satellites or minor com- 
panions, can be considered part of the halo of the large host 
galaxy. 

In the following subsections, we discuss only the ten galax- 
ies selected for this study (see column 8 of Table ID . Detaile d 
images of the galaxies are presented in K acprzak et al.l (l2007h . 
A summary of the impact parameters, galaxy maximum rota- 
tion velocities, and GIM2D model inclinations and position 
angles (PA = angle between the galaxy major axis and the 
quasar line of sight) is listed in Tabled Galaxy redshifts will 
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Fig. 2.— (a) A 30" x 20" F702W WFPC-2/HST image of the quasar field Q0002+051.The ESI/Keck sHt is superimposed on tlie image. Tlie '+' and '-' on 
the slit indicate the positive and negative arcseconds where 0" is defined at the target galaxy center. The targeted galaxy is a compact galaxy with z = 0.8518 
at impact parameter D = 25.9 kpc. Two other galaxies have been identified in this field and their redshifts are indicated. — (b) The z = 0.8518 galaxy rotation 
curve and the HIRES/Keck absorption profiles aligned with the galaxy systemic velocity. — (c) Same as (a) except for the Q0229+ 131 quasar field. The targeted 
galaxy is a low inclination spiral with z = 0.4167 at D = 37.5 kpc. — (d) Same as (b) except the z = 0.4167 galaxy in the Q0229 + 131 field (Mgl A2853 absoiption 
is not shown because it is blended with SilV A1394 associated with a z = 1.9024 CiV absorber). Note that the IVIgll absoiption resides fully to one side of the 
galaxy systemic velocity and also aligns with one arm of the rotation curve. 

pears Gl is not predominantly rotating. The projected veloc- 
ity shear is 49 km s"'. Displayed in the lower panels are the 
Mgll and Mgl absorption profiles, shown on the same veloc- 
ity scale. The mean optical depth absorption redshift is offset 
by -66 km s"' from the galaxy systemic velocity. This Mgll 
absorption system is the most kinematically complex absorp- 
tion profile in our sample. Four absorption sub-systems are 
spread out over a total velocity of 475 km s"'. The blue- 
shifted extreme wing of Gl's velocity shear is partially consis- 
tent with the velocities of the dominant saturated component 
of the Mg II. However, it is inconsistent with the Mg II gas at 
large positive velocities. The Mgl gas, which may trace more 
neutral and higher column density gas, has velocities is also 
inconsistent with the velocities of Gl. 



only be quoted to four significant figures from here on for sim- 
pUcity. We will later discuss kinematic halo models in § |4] 

3.1. Q0002+051 Gl 

A WFPC-2/HST image of the Q0002 + 051 field is shown 
is FigurelJJ?. The z = 0.8518 g alaxy, Gl, was targeted for this 
study. It was first assumed by ISteidel. Dickinson. & Perssoa 
(1 19941) to be the absorbing galaxy associated with the Mgll 
absorption at z = 0.851407 given its colors and its proximity 
to the quasar. We report the first redshift confirmation of this 
galaxy. The galaxy redshift was identified by an [Olll] and 
a weak [Oil] emission line (the [Oil] line was used only for 
the determination of the galaxy redshift and provided no kine- 
matic information). In Figure |2j5(, we note that the ESI/Keck 
slit was placed across both Gl and G2 (where '+' and '-' on 
the slit indicate the positive and negative arcseconds relative 
to the center of the galaxy, respectively). We did not detect 
any identifiable emission lines from the galaxy G2. 



The Mg II absorption was first reported by Bechto ld et al 
and then confirmed 

The HIRES/Keck spectrum 



(11984) 
1992) 



by ISteidel & Sargent 



was originally 

presented by IChurchill. Steidel. & Vogll d 19961) and has 
W,(2796)= 1.119 A. 

The Gl galaxy has a compact morphology. It has an ab- 
solute Z?-band magnitude of Mb = -21.3 which translates to 
a luminosity Lb = 0.92L|. The galaxy has inclination an- 
gle i = 38° and PA = 43°. Gl is at a projected distance of 
D = 25.9 kpc from the quasar line of sight. In Figure [JJ?, the 
spatial radial velocity, as derived from [Olll], is shown. It ap- 



3.2. Q0229+131 Gl 

A WFPC-2/HST image of the Q0229+ 131 field is shown 
in Figure lit. The z = 0.41 67 galaxy, Gl, was i nitially spec- 
troscopically identified by [Bergeron & Boissl (1199 lb . The 
first detection o f Mgll a bsorptio n at the galaxy red shift 
was reported by ISargent. Boksenbe rg, & Steidell (I l988 l). A 
HI RES/Keck spectrum of this ab sorber was first presented 
bv lChurchill. Kacprzak. & Steidell f2005i) and has W,(2796) = 
0.816 A. 

This moderately inclined, / = 58°, bulge dominated galaxy 
has two large grand design spiral arms similar to a local Sb 
galaxy. One of the arms contains a bright Hll region. The 
galaxy position angle is PA = 22°. Gl has an absolute B-band 
magnitude of Mb = -20.8 which translates to Lb = 0.74L^. Gl 
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Fig. 3. — (a) Same as Figure|2j( except for the quasar field Q0450- 132. The target spiral galaxies Gl and G2 are possibly interacting galaxies, indicated by 
the single-sided tidal tails. Gl is at a redshift of z = 0.4942 at D = 50.1 kpc. GI is at a redshift of z = 0.4931 at D = 62.7 kpc. — (b) Same Figure|2f) except Gl 
and G2 in the Q0450— 132 field. — (c) Same as (a) except the z. = 0.4838 galaxy, Gl, in the Q0454— 220 field. The targeted galaxy is a spiral galaxy at an impact 
parameter D = 107.9 kpc. — (d) Same as (b) except for the Gl galaxy in the Q0454 — 220 field. 



is at a projected distance of D = 37.5 kpc from the quasar line 
of sight. 

In Figure |2]i, the rotation curve, derived from [Oil], H/3, 
and [Olll] emission lines, is shown. The asymmetric rota- 
tion curve has a maximum observed velocity of 280 km s~' . 
The asymmetry may be a result of the asymmetric spiral 
arms. Below the rotation curve, the Mgll AA2796,2803 ab- 
sorption profiles are shown on the same velocity scale. The 
mean absorption redshift is offset by +135 km s~' from the 
galaxy systemic velocity. The Mgl profile is not shown 
here since it is completely bended with a four compo- 
nent Si IV A 1394 complex froni a z = 1.9024 CiV absorber 
(ISargent. Boksenberg. & Steidell fT988h . This Mgll absorp- 
tion system has a velocity spread of roughly 1 12 km s"' with 
a single cloud residing roughly +280 km s"' from the main 
component. The main absorption component aligns with the 
redshifted wing of the galaxy rotation curve arising from the 
spiral arm nearest the quasar The galaxy rotation velocities 
are not consistent with the outlying cloud. 

3.3. Q0450-132G1,G2 

A WFPC-2/HST image of the Q0450- 132 field is 
shown in Figure Two galaxies were placed along 

the same slit. The redshifts of this double pa i r wer e 
initially obtained by ISteidel. Dickinson. & PerssonI (|1994|) . 
This absorption system was discovere d during the survey of 
ISargent. B oksenberg. & S teidel ( 1988) and the HIRES/Keck 
spectrum was first presented by Churchill ( 1997). The Mgll 
absorption has an equivalent width of Wr(2796) = 0.674 A. 

The galaxy closer to the quasar, Gl, has a redshift of z = 
0.4941 and is at an impact parameter of D = 50.1 kpc. The 



galaxy has an inclination of i = 66° with a PA= 42°. Gl has 
an absolute magnitude of Mb = -19.7 which translates to a 
Lb = Q.25Lg galaxy. The galaxy further from the quasar, G2, 
has a redshift of z = 0.4931 and has an impact parameter of 
62.7 kpc. The galaxy has an inclination of / = 75° with a 
PA = 54°. G2 has an absolute magnitude of Mb = -19.7 which 
translates to Lb = 0.26Lg. The galaxies are spatially separated 
by only 14 kpc projected, and by a line of sight velocity of 
201 km s"'. Both galaxies appear to have single-sided tidal 
tails and show strong morphological evidence of a previous 
interaction or harassment. 

The rotation curves of both galaxies, as derived from [Oil], 
H7, up, [Olll], Ha, and [Nil], are presented in Figure [S?. 
The maximum observed rotation velocities for Gl and G2 are 
98 km s"' and 47 kms"', respectively. Below the rotation 
curve, the Mg II and Mg I absorption profiles are shown on the 
same velocity scale. The mean absorption redshift is offset 
by -33 km s"' from the galaxy systemic velocity of Gl. The 
Mg II absorption is a single kinematic region with a velocity 
width of 103 km s"'. The absorption aligns, in velocity, with 
the blueshifted wing of the rotation curve of Gl. The Mgl 
aligns more closely to the systemic velocity of Gl . The galaxy 
G2 does not have much rotation since the slit position angle 
is close to its minor axis. 

3.4. Q0454-220G1 

A WFPC-2/HST image of the Q0454-220 field is shown 
in Figure [3]:. The z = 0.4838 galaxy, Gl, was targeted 
for this study. Gl was spectroscopically identified by 
|Chen et al.' ('1998|). The M gll absorption was first repo rted by 
iBergeron & KunM d 19841) and then confirmed by Tvtler et alj 



Kacprzak et al. 



Q0836+113 zgal=0.7868 



Q0836+113 zgal=0.7868 



I [oil] A:i7a7 



z(Gl)=0.7868-l 



-200 




D(G1)=27 kpc 

m 

QSO 

■• — z(G2)^0.4829 



Q1127-145G3 zgal-0.3285 



z(GiM.3132 



QSO||f 



z(G2)^0.3124 



4 



z(G3)=0.3285 ^ A 
D(G3)=92 kpc ^ 



* <- z(G4)^0.312l 



Ql 127-145G3 zgal=0.3285 



fNII] A65B5 
Ha X6^)B3 



2 - 



m 



-200 

1 .5 \^ 



\---- 
-200 



-100 100 200 

Mgll X2796 - 



-100 100 200 



Mgl A2053 



-100 

vel (km/a) 



100 200 



Fig. 4. — (a) Same as Figure|2j( except for the quasar field Q08364- 1 13. The target spiral galaxy Gl is at a redshift of z = 0.7868 at Z) = 26.9 kpc. — (b) Same 
Figure|2j) except Gl in the Q08364- 1 13 field. — (c) Same as (a) except the z = 0.3285 galaxy, Gl, in the Ql 127- 145 field. The targeted galaxy is a spiral galaxy 
at an impact parameter D = 91.4 kpc. — (d) Same as (b) except for the Gl galaxy in the Ql 127 — 145 field. 

sided tidal tails. This almost edge-on galaxy has an incli- 
nation of = 78° and M = 57°. Gl has an absolute B-band 
magnitude of Mb = -20.5 which translates to Lb = 0.39Lg. 
The galaxy is at a projected distance of D = 26.9 kpc from the 
quasar 

The spatial radial velocity of Gl, as derived from [Oil], 
is presented in Figure |4}7. The data suggest that Gl exhibits 
more of a global shear than rotation. The maximum observed 
shear velocity is ^ 42 km s~' . Below the velocity curve, 
the Mg II and the Mg I absorption profiles are shown on the 
same velocity scale. The mean absorption redshift is offset 
by -16 km s"' from the galaxy systemic velocity. The Mgll 
absorption is a single kinematic component having a velocity 
spread of roughly 282 km s"' . Most of the component is com- 
posed of highly saturated clouds. From the Mgl profile, one 
can resolve the individual clouds that are saturated in Mgll. 
The galaxy velocities are consistent with the velocities of the 
absorbing gas. However, there is a large amount of Mgll gas 
that has greater velocities than those of the galaxy. The bulk 
of the more neutral/high column density gas, as indicated by 
Mgl, is at the galaxy systemic velocity. 



d Hig). The HIRES/Keck sp ectrum was originally presented 
by IChurchill & Vogti ( 120011) . The Mgll equivalent width is 
W^,(2796) = 0.426 A. 

Gl is a spiral galaxy that has a perturbed morphology with 
one extended spiral arm. It has a compact bulge and several 
bright Hll regions. The galaxy has an inclination of / = 41° 
and a position angle of PA = 16°. It has an absolute B-band 
magnitude of Mb = -21 .0 which ti-anslates to Lb = 0.90L^. Gl 
is at a projected distance of D = 107.9 kpc from the quasar line 
of sight. 

The rotation curve, as derived from [Oil], H/3, [Olll], Ha, 
and [Nil], is presented in Figure [3]i. The rotation curve flat- 
tens out at a maximum observed velocity of 138 km s~'. Be- 
low the rotation curve, the Mg II and the Mg I absorption pro- 
files are shown on the same velocity scale. The mean ab- 
sorption redshift is offset by -98 km s"' from the galaxy sys- 
temic velocity. This Mgll absorption is a single kinematic 
system having a velocity spread of roughly 93 km s"'. The 
blue shifted component of the rotation curve has velocities 
that are consistent with the Mg II and Mg I absorption profiles. 



3.5. Q0836+113G1 

A WFPC-2/HST image of the Q0836+ 113 field is shown 
in Figure lU;. The z = 0.7 868 galaxy, Gl, was init ially spec- 
troscopically identified by iLowenthal et akl (Il995h . The first 
M gll absorption d etection at the galaxy redshift was reported 
by iTurnshek et al.1 ([1989). The HIRES/Keck absorption pro- 
files of this absorber is first presented in this work (FigurelH?) 
and has W,(2796) = 2.15 ±0.02 A. 

The galaxy Gl appears to be an edge-on spiral with an 
asymmetric brightness profile. The morphology is similar to 
that of the two galaxies in Q0450- 132, which have single- 



3.6. Ql 127-145 G3 

A WFPC-2/HST image of the Ql 127- 145 field is shown 
in Figure HJ:. We present a newly identified galaxy, G3, at 
z = 0.3285. The galaxy was identified by Ha and [Nil] emis- 
sion lines. The absorption was also recently discovered, and is 
presented here in FigureHlf. This weak system has an equiv- 
alent width of W,(2796) = 0.029 A. 

The face-on galaxy has an inclination of / = 1° with a 
PA ^ 69° (the PA is highly uncertain since / ~ 0). The galaxy 
has a large bar with a sizable bright bulge; similar to a local 
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Fig. 5.— (a) Same as Figure |2j; except that this is a 30" X 30" F814W WFPC-2///5r image of the quasar field Ql 127- 145. The two targeted spiral galaxies 
Gl and G2 have redshifts of z = 0.3132 and z = 3124, respectively. Gl and G2 are at impact parameters of D = 45.6 kpc and D = 81.0 kpc, respectively. — (b) 
Same as Figure|2j) except Gl in the Ql 127—145 field. — (b) Same as Figure|2j) except G2 in the Q1127 — 145 field. 



SBb galaxy. It has absolute magnitude of Mb = -20.51 which 
translates to Lb = 0.60L^. The galaxy is at a projected distance 
of D = 91 .4 kpc from the quasar 

The rotation curve of G3 is presented in Figure HJ/. The 
maximum observed rotational velocity is ^ 80 km s"' . Given 
that significant rotation is observed, the galaxy is most likely 
not completely face-on as the GIM2D model inclination sug- 
gests. Below the rotation curve, the Mgll and the Mgl ab- 
sorption profiles are shown on the same velocity scale. The 
absorption redshift is offset by -46 km s"' from the galaxy 
systemic velocity. The Mg II absorption contains two separate 
single cloud components. Both clouds have a velocity spread 
of 15 km s~' and are separated by 125 kms"'. No sig- 
nificant Mgl (3(7, Wr < 0.003 A) is detected. The projected 
galaxy rotation velocities are consistent with both cloud ve- 
locities; each cloud aligns with each side of the rotation curve. 

3.7. Q1127-145G1,G2 

The history of the Q 11 27- 145 z = 0.3 12710 absorption sys- 
tem is q uite complex. The Mgll ab sorption was initially de- 
tected by'Beraeron & Boisse' ('1991') and was determined to be 
a DLA (see Rao & Turnshek 2000) since H ST UV data show 
a damped Lya with Nhi = 5.1 x 10^' cm"^ (ILane et al.lll998h . 
The equivalent width of the system is Wr(2796) = 1 .773 A. 

The true identity of the absorbing galaxy has been a topic 
of debate in the literature. A WFPC-2 /HST image of the 
Ql 127 - 145 field is shown in Figure |5]j!. iBergeron & Boissa 
(1199 1 ) spectroscopically identified Gl and G2 to be at the 
redshift of the absorption. Gl exhibits strong emission lines 
and G2 has no detectable emission lines. Gl was assumed 
to be the absorbing galaxy, since it is closer to the quasar 
and has significant star formation (.Bergeron & Boisse 1991.) . 



iLane et all (1 19981) later spectroscopically identified G4 via an 
[Olll] doublet with z = 0.3121, which was also consistent with 
the absorption redshift and was assu med to be the absorber 
only due to its proximity to the quasar. iLane e t al. ( 1998) state 
that it is also possible that the three galaxies may have un- 
dergone a strong interaction where the absorption could arise 
from tidal debris. In Figure |5]2, it is apparent from the tidal 
disturbances, that Gl has undergone interactions in the past. 
G2 also exhibits some tidal material to the north of the galaxy 
( which is les s appar ent in the figure). 

iRao et al.1 (l2003h and iNestor et al.l (l2002h suggest, from 
ground-based multi-band imaging, that the low surface 
brightness emission detected around the quasar (see Fig- 
ure ISjj) could arise from a foreground low surface bright- 
ness galaxy at the absorption redshift. However, it is pos- 
sible that the low surface brightness signal is coming from 
the quasar host galaxy at Zem = 1-18. The background quasar 
is radio-loud and has strong X-ray emission. An X-ray 
j et extends north east 236 kpc f rom the quasar, projected 
dSiemiginowska et al.ll2002l |2007|) . There is a diffuse halo, 
both detected in radio and X-ray, around the host quasar It 
is not yet clear whether the X-ray halo is real or a result of 
blurring from the instrument PSF (Siemiginowska, A. 2007, 
private commutation). Th us, it is possible that the low sur- 
face brightness detected by Rao et al. (2003) may be from the 
background quasar. (Chun et al. 2006) found a possible un- 
derlying galaxy 0.6" from the quasar. Again, this may be a 
foreground galaxy, or structure from the quasar host galaxy 
which is commonly observed (e.g., Bahcall et al. 1997). 

Galaxies Gl and G2 were targeted for this study. These 
galaxies are spatially separated by only 55.3 kpc projected, 
and by a line of sight velocity of 183 km s"' . Gl is an almost 
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Fig. 6. — (a) Same as Figure|2}; except for the quasar field Q2206- 199. 
same FigureiJ) except Gl in the Q2206- 199 field. 

edge-on (/ = 82°) spiral that displays asymmetries on both 
sides of the galaxy. The galaxy has PA = 21°. It has absolute 
magnitude of Mb = -20.40, which translates to Lb = 0.54L^. 
The galaxy is at a projected distance of D = 45.6 kpc from the 
quasar line of sight. G2 is an interesting galaxy; it has a major 
dust lane and a large bulge. Given that we detected no emis- 
sion lines, this galaxy could either be classified as Sa or as an 
early-type SO galaxy. G2 also has a tidal disturbance along 
the major axis of the galaxy towards the north. The galaxy 
has an inclination of / = 73° and PA = 61°. It has absolute 
magnitude of Mb = -20.57 which translates to Lb = 0.63L^. 
The galaxy is at a projected distance of D = 8 1 .0 kpc from the 
quasar. 

The rotation curve of Gl, obtained from the [Oil], Ha, and 
the [Nil] doublet, is presented in Figure 13?. The maximum 
observed rotational velocity is 204 km s"' . The rotation curve 
of G2, obtained from the Nal absorption doublet, is presented 
in Figure |5];. The maximum observed rotational velocity of 
G2 is 90 km s"'. Below both rotation curves, the Mgll and 
the Mg I absorption profiles are shown on the velocity scale 
defined by Gl. The mean absorption redshift is offset by 
-112 km s"' from the systemic velocity of Gl. The Mgll 
absorption can be broken up into two kinematic components. 
The first large saturated component has a velocity spread of 
roughly 235 kms"'. The second component, blue-ward of 
the main component, contains a few weak clouds and has a 
velocity spread of roughly 68 km s~' . In the Mg I profile, one 
can resolve the individual clouds that are saturated in Mgll. 
Almost all of the Mg I gas is aligned with the saturated Mg II 
component. Only a very weak Mgl cloud is detected in a 
second kinematic component. The main component of the 
Mgll gas aligns with the blue-ward wing of the Gl rotation 
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target spiral galaxy Gl is at a redshift of z= 1.10166 at Z) = 104.6 kpc. — (b) 

curve, which is on the side toward the quasar. The Mgll ab- 
sorption velocities start at the systemic velocity of Gl and is 
centered on the maximum galaxy rotation velocity. Since the 
rotation curve of G2 is derived from Nal absorption lines, it 
is likely that the maximum rotation speed of the galaxy ex- 
tends to larger velocities then detected. The velocities of G2 
are consistent a portion of the large saturated component of 
the Mgll and the small blue-ward clouds. 

3.8. Q2206-199G1 

A V^FPC-2/HST image of the Q2206- 199 field is shown 
is Figure The z = 1.0166 galaxy, Gl, was targeted 

for this study. The galaxy w as spectroscopically identi- 
fied by iBergeron. Cristiani. & S haver ( 1992). The Mgll ab- 
sorption was first reported bv lSargent. Boksenberg. & Steidell 
(198 8). The HIRES/Keck spe ctrum was originally presented 
by lProchaska & Wolf d (Il997h . The Mgll equivalent width is 
Wr(2796)= 1.057 A. 

The spiral galaxy Gl is quite unusual in brightness and mor- 
phology. The galaxy has an absolute magnitude of Mb = -23.8 
which translates to Lb = 8.0L^. Galaxies with super-L* lu- 
minosities are quite rare and represent only a few percent of 
the galaxy population. The galaxy spiral structure is tightly 
wound with a large leading arm. There appears to be numer- 
ous Hll regions. The bulge is compact and offset from the 
isophotal center Note that at the galaxy redshift, the rest- 
frame mean wavelength of the F702W filter is around 3440 A, 
roughly rest-frame U-band. The galaxy inclination is / = 57° 
with a PA = 67°. The galaxy is at a projected distance of 
D = 104.6 kpc from the quasar. 

The rotation curve, obtained from the [Oil] doublet, is pre- 
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sented in FigurelSI?. The maximum observed rotational veloc- 
ity is 40 km s~' . Given the low observed rotation velocities, 
the galaxy is likely more face-on than the G1M2D model in- 
clination suggests. The galaxy morphology is asymmetric and 
unusual which makes it difficult to determine the inclination. 
Below the rotation curve, the Mgll and the Mgl absorption 
profiles are shown on the same velocity scale. The mean ab- 
sorption redshift is offset by 73 km s~' from the galaxy sys- 
temic velocity. The Mgll absorption can be broken up into 
two kinematic components. The first one has velocity spread 
of roughly 150 km s"' and is mostly saturated. The second 
component is a single cloud with a velocity spread of roughly 
28 km s"' that is offset 198 km s"' red-ward of the main com- 
ponent velocity center From the Mg I profile, one can resolve 
the individual clouds that are saturated in Mgll. Most of the 
Mgl gas is shifted 100 km s"' from the systemic velocity of 
the galaxy. The red-ward arm the the rotation curve aligns 
with only a small portion of velocity of the main Mg II com- 
ponent. However, there is a large amount of Mg II gas that has 
greater velocities than those of the galaxy. The bulk of the 
more neutral/high column density gas, as indicated by Mgl, 
has velocities greater than that of the observed galaxy rota- 
tional velocities. 

3.9. Summary I: Observational Kinematic Comparisons 

In our direct comparison of galaxy disk and halo gas kine- 
matics, traced by Mg II and Mg I absorption, we find the fol- 
lowing: (1) in all ten cases, the observed galaxy rotation ve- 
locities show substantial overlap with the bulk of the absorp- 
tion velocities. (2) in seven of ten cases, the Mgll and Mgl 
absorption resides to one side of the galaxy systemic velocity. 
In the remaining cases, (Q0002 + 051 Gl, Q0836+113 Gl, 
and Ql 127 - 145 G3) absorption resides on both sides of the 
galaxy systemic ve locity. Our findings are similar to those of 
ISteidel et al.l(l2002h . even though their sample targeted highly 
inclined and edge-on disk galaxies with PA ^ 0. Here we 
have attempted to select galaxies with a range of inclination 
and position angle with respect to the quasar line of sight (see 
Table©. 

To see if there are differences in the Mg II absorption pro- 
files as a function of inclination, we have separated the galax- 
ies into two inclination bins with / < 60° (five g alaxies) and 
i > 60° (six galaxies). The galaxies from the Steidel et~an 
([2002) sample are included here. For this comparison, mul- 
tiple galaxies that can be associated with a single absorption 
system (such as Gl and G2 in Q0450- 132 and 1127-145) 
were removed since we cannot confidently know whether one 
or both galaxies host the Mg II absorption. 

In Figure |2l we show the co-added Mg II absorption pro- 
files. The data are plotted as an absolute velocity difference 
from the galaxy systemic velocity. As seen in Figure [TJj, the 
combined spectra of the galaxies with / < 60° shows that ab- 
sorption resides between 0-180 km s"' with a peak in the 
optical depth at v ^ 100 km s"'. In Figure iTf), the combined 
spectra show that galaxies with / > 60° are associated with 
Mgll absorption with smoothly varying optical depth spread 
over ~ 300 km s"' . A Kolmogorov-Smirnov (K-S) test re- 
veals the probability of the two optical depth distributions be- 
ing drawn from the same sample is = 0.000873. This 
rules out the null hypothesis of similar Mgll optical depth 
distributions as a function of inclination at the 99.91% con- 
fidence level (3.3cr). The data are suggestive that Mgll ab- 
sorption velocity spread may be a function of galaxy orien- 
tation; galaxies with higher inclination have Mg II absorption 
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Fig. 7. — The combined Mgll absorption spectra for galaxies separated 
into two inclination bins. The spectra are plotted versus absolute velocity dif- 
ference from the galaxy systemic velocity. The five galaxies of Steidel et al] 
(2002) are included. The two galaxy pairs, Gl and G2 of Q0450- 132 and 
Ql 127 - 145, are excluded here. — (a) The summed spectra for five galaxies 
with i < 60° . — (b) The summed specfi'a for six galaxies with i > 60° . 

with larger velocity spreads and more evenly distributed opti- 
cal depths (on average at any given velocity). However, since 
the number of galaxies per inclination bin is small, we need 
to acquire a larger sample to see if the trend holds. 

A trend with inclination might be expected if the absorbing 
gas kinematics is well represented by a monolithic rotating 
halo. Mgll ab sorption velocity spr eads in four out of the five 
galaxies in the lSteidel et al.l (12002*) sample were shown to be 
consistent with a monolithic rotating halo model that allowed 
for lagging rotation with increasing height above the disk. 
Given the apparent kinematic trend with inclination, we in- 
vestigate whether the lagging halo model (ISteidel et alj|2002l) 
can successfully predict the Mg II absorption velocity spreads 
of our ten galaxies. 

4. GALAXY KINEMATICS AND HALO-DISK MODELS 



We apply the simple halo model of ISteidel et alj (l2002h to 
our systems in order to determine whether an extended disk- 
like rotating halo is able to reproduce all or most of the ob- 
served Mg II absorption velocity spread. The model is a co- 
rotating disk with velocity decreasing as a function of scale 
height. 

The line of sight velocity, v/,„, predicted by this disk halo 
model is a function of the measurable quantities D, i, PA (the 
angle between the galaxy major axis and the quasar line of 
sight), and v^ax, which is the maximum projected galaxy ro- 
tation velocity. 
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D sin PA 

yo = and p = DcosPA, 

cos/ 

where the free parameter, /i,,, is the lagging gas velocity scale 
height and where y is the projected line of sight position above 
the disk plane. The parameter represents the position at 
the projected mid-plane of the disk. The range of y values is 
constrained by the model disk-halo thickness, //eff, such that 
yo- Hsfftani <y <yo +//efftan/. The distance along the line 
of sight relative to the point were it intersects the projection 
of the disk mid-plane is then Di,„ = (y-yo)/ sin/. There are 
no assumptions about the spatial density distribution of Mg II 
absorbing gas, except that H^ff is the effective thickness of the 
gas layer capable of giving rise to absorption. 

In order to maximize the rotational velocity predicted by the 
model, we assume /;,, = 1000 kpc, which effectively removes 
the lagging halo velocity component (such that the exponen- 
tial in Equation[T]is roughly equal to unity). 

In Figure [8]3-/i, we show the Mgll absorption profiles for 
each galaxy, where the shaded regions indicate detected ab- 
sorption. Below each absorption profile is the disk halo model 
velocities as a function of Dios derived for each galaxy (solid 
line) using Equation[T]and parameters in Table |6] Recall that, 
at Djgs = kpc, the model line of sight intersects the projected 
mid-plane of the galaxy. The dashed curves represent the 
disk halo model velocities derived from the combination of 
the minimum and maximum uncertainties in the PA and /. In 
some cases (see Figure |8l7) the values of the PA and / are well 
determined such that the dashed curves lie on the solid curves. 
The model also predicts the line of sight position, D/oj, of the 
halo gas at each velocity, v/g,. 

The disk halo model is successful at predicting the ob- 
served Mgll absorption velocity distribution when the solid 
(or dashed) curves span the same velocity spread as that of 
the Mgll absorption gas. The model curves must occupy the 
full shaded region to be 100% successful. If this is not the 
case, one can conclude that disk-like halo rotation is not the 
only dynamic mechanism responsible for the Mg II kinemat- 
ics. In the following subsections we discuss the disk model of 
the individual galaxies. 

4.1. Q0002+051 Gl 

The galaxy Gl exhibits a low level velocity shear. Given 
the velocity spread of the gas (^ 475 km s"'), it is impossi- 
ble for the bulk of the absorption gas to be consistent with 
the observed velocities of Gl. In Figure [8]3, we see that the 
galaxy disk halo model is counter rotating with respect to the 
dominate saturated Mg II component. There is no overlap be- 
tween the predicted halo model velocities with those of the 
Mg II and Mg I absorption. Even if the galaxy had a highly 
significant velocity shear, the bulk of the Mg II clouds would 
not be consistent in velocity space. Given the number of high 
velocity components, it is unclear that this absorption profile 
represents a gravitationally bound gaseous galactic halo. 

4.2. Q0229+131 Gl 

The galaxy Gl has an asymmetric rotation curve with the 
largest rotation velocity observed in the direction of the quasar 
line of sight. In Figure ISj?, we see that the disk halo model 
velocities are consistent with the bulk of the Mg II absorption. 
Thus, extended disk-like halo rotation could be invoked to ex- 
plain most of the observed halo gas velocities. However, there 



TABLE 6 
Galaxy Disk Model Input Values 



QSO Field 


Gtilflxy 
ID 


D 
(kpc) 


(km/s) 


(deg.) 


PA 
(deg.) 


Q0002 + 051 


Gl 


25.9 ±0.5 


-49 




43!" 


Q0229+13I 


Gl 


37.5 ±0.5 


-281 


58!f 


22^2 


Q0450-I32 


Gl 


50.1 ±0.4 


-98 


66;3 


42+2 


Q0450-I32 


G2 


62.7 ±0.7 


47 


75;2 


54+2 


Q0454 - 22() 


Gl 


107.9 ±0.8 
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41!^ 




Q0836+I13 


Gl 


26.9 ±0.9 


-42 


78!] 




QI127-I45 


Gl 


45.6 ±0.3 


204 


82:« 


,1+0 


QI127-I45 


G2 


81.0±0.3 


-90 


73!^ 


61^2 


QI127-I45 


G3 


91.4±0.2 


-80 




69!?^ 


Q2206-I99 


Gl 


104.6 ± 1.4 


-40 




61% 



remains a small single Mgll halo cloud ^ 400 km s"' from 
systemic that that cannot be explained by halo rotation alone. 
This suggests other dynamic processes give rise to some of 
the Mg II absorption. 

4.3. Q0450-132 Gl, G2 

Galaxies Gl and G2 are potentially interacting galaxies, as 
evident from their morphologies and strong emission lines. 
It is possible that these interactions are an effective mech- 
anism in producing extended Mgl l abs orption in th e halo 
jBowen et al.l l 1995; Kacp rzak et alj|2007l: iRubin et al.l l2009^. 
In Figure [SJ:, we plot the disk halo models for both galax- 
ies. Gl, the galaxy closest to the quasar, has model halo gas 
kinematics that are counter-rotating with respect to the Mgll 
absorption. G2, on the other hand, has modeled halo veloci- 
ties that are consistent with those of the Mg II absorption. The 
halo model of G2 is also consistent with the bulk of the Mgl. 
Given that Gl was observed along the major axis and G2 was 
not, if we assume that G2 had comparable rotation speeds as 
Gl, the halo model velocities would overlap with most of the 
absorption velocities. We will discuss the difficulties of dis- 
entangling these multiple galaxy systems in § 14.91 

4.4. Q0454-220G1 

Galaxy Gl has a symmetric rotation curve that completely 
flattens out at the maximum velocity of 138 km s"'. In Fig- 
ure[8]i, we see that the disk halo model has velocities consis- 
tent with the bulk of the Mgll absorption velocities. They are 
also consistent with the Mgl absorption velocities. However, 
there is an inconsistency of ^ 40 km s"' between the halo 
model and Mgll absorption velocities. Thus, the halo model 
is unable to reproduce the total observed spread of Mg II ab- 
sorption velocities. 

4.5. Q0836+113G1 

The galaxy Gl exhibits minimal rotation; the velocities are 
more indicative of a global shear The galaxy systemic ve- 
locity is centered roughly in the middle of the Mg II and the 
bulk of the Mgl absorption profiles. In Figure[8^, we see that 
the halo model also shows little rotation, roughly 50 km s"' . 
If the model was a true representation of the halo, then more 
than 50% of the absorbing gas has velocities inconsistent with 
disk rotation that are larger than the model velocities. Even if 
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Fig. 8. — The Mgll absorption profiles and the disk model velocities as a function of D/oj (solid curve) are shown for each galaxy in the top and bottom 
panels, respectively. The Mgll absoiption velocities are shaded in. The solid curve is computed using Equation [T] and the values from Table|6] The dashed 
curves are models computed for the maximum and minimum predicted model velocities given the uncertainties of and PA. The disk model is successful and 
reproducing the observed absorption velocities in the solid curve overlaps with the entire shaded region. D;„j is equal zero when the quasai' line of sight intersects 
the projected mid-plane of the galaxy. The panels ai'e as follows; (a) Q0002 + 051 Gl, Gl (b) Q0229+ 131 Gl, (c) Q0450- 132 Gl and G2, (d) Q0454-220 Gl, 
(e) Q0836+ 1 13 Gl, (f) Ql 127-145 G3, (g) Ql 127- 145 Gl and G2, and (h) Q2206- 199 Gl. 
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Gl has a more significant velocity shear, the model would 
still not be able to explain the gas blue-ward of the galaxy 
systemic velocity. Here, the models fails to predict the bulk 
of the absorption velocities. 

4.6. Q1127-145G3 

The galaxy G3 appears face-on, however it exhibits a max- 
imum rotation of 80 km s"' . The model inclination of 1 ° 
is likely incorrect given the observed rotation velocities. In 
Figure [8]f , the disk halo model for G3 exhibits little line of 
sight velocity 20-50 km s"') given its orientation with re- 
spect to the quasar line of sight. Given the model parameters 
listed in Table |6] we varied the galaxy inclination, such that 
< / < 45°, in an attempt to reproduce the observed absorp- 
tion velocities. Even with / = 45°, the halo model fails to re- 
produce the observed Mg II absorption velocities. For the disk 
halo scenario, the halo gas is expected to reside to one side of 
the rotation curve. The nature of this Mg II absorption profile 
is interesting; two weak clouds separated by 125 km s"' . It is 
likely that these clouds could arise in either a patchy diffuse 
halo or the line of sight is intercepting small scale structure 
near the galaxy halo. In any case, the model is unable to re- 
produce the observed Mgll absorption velocities. 

4.7. Ql 127-145 G1,G2 

Galaxies Gl and G2 are potentially interacting galaxies, as 
evident from their morphologies. These galaxies appear to be 
in a small group including G4. In Figure [8^, we plot the disk 
halo models for Gl and G2. The model for Gl has velocities 
that are consistent with up to 200 km s"' spread of the satu- 
rated component of the Mg II and all of the Mg I absorption. 
The disk halo model of G2 is counter-rotating with respect to 
Gl as viewed from the quasar line of sight. The model veloc- 
ities are consistent with the remaining absorption of the sat- 
urated component which was not covered by Gl. The model 
velocities of both Gl and G2 overlap ^ 75 km s"'. The sec- 
ond weaker kinematic component, ^ 300 km s"' blue-ward 
of systemic velocity of Gl, cannot be explained given the pre- 
dicted halo velocities. It is possible that the saturated com- 
ponent of the Mgll absorption could arise from either Gl, 
G2, or both. The weaker component may arise from tidal 
debris stirred up by the apparent interactions. It is also pos- 
sible that some of the Mg II absorption is associated with G4. 
It is likely that the saturated component is associated with a 
nearby galaxy since it is a DLA and many other low ions have 
also been detected. We will discuss the difficulties of disen- 
tangling these multiple galaxy systems in § 14.91 In any case, 
the model is mostly successful, except that it does not repro- 
duce all of the Mgll absorption velocities. 

4.8. Q2206-199G1 

The Gl galaxy is the largest in our sample, and the second 
furthest away from the quasar line of sight. In Figure [8^, the 
disk halo model shown has very little line of sight velocity 
40 km s"'). It is clear that the model of this moderately in- 
clined galaxy does not reproduce the observed absorption ve- 
locities. In fact, the bulk of the Mgll and Mgl has velocities 
^ 50 km s"^ greater than the galaxies maximum observed ro- 
tation velocity. This is peculiar, since the dominant saturated 
component is commonly expected to be associated with the 
galaxy disk, yet the kinematics here suggest otherwise. This 
galaxy-absorber pair is a another example demonstrating that 
disk-like halo rotation cannot be the only mechanism driving 
the kinematics of halo gas. 



4.9. Mgll Absorption from Galaxy Pairs/Groups 

Since Mgll absorbers were first associated 
with galaxy halos, it has been common practice 
to associate one galaxy with an absorption sys- 
tem at a given redshift (e.g., Bergeron & Boissa 

19911: ISteidel. Dickinson. & Persso n" Il994t 

Churchill. Steidel. & Vogt 1996; Guillemin & BergeronI 
1993). It is now becoming evident th at Mgll absorption 
also arises in small groups of galaxies dBowen et al.l 119951; 
Churchill & Charlton 1999) and even clusters jLopez et^ 
2008). In this paper, we present two such examples: 
Q0450-132 Gl, G2 and Q1127-145 Gl, G2, G4. The 
pair of galaxies in Q0450- 132 are close to each other in 
both projected distance and velocity. Both galaxies have 
morphological evidence (one sided tidal tails) that is clas- 
sically associated with interacting/merging galaxies. The 
Ql 127- 145 field contains three galaxies at similar redshift 
(we have recently spectroscopically identified two additional 
galaxies with similar redshifts within D = 250 kpc from 
the quasar [Kacprzak et al., in prep]). The two galaxies 
studied here, Gl and G2, have evidence of morphological 
perturbations and extended tidal material. 

The Q0450- 132 and Ql 127- 145 galaxy pairs, have ob- 
served rotation velocities that overlap with those of the Mg II 
absorption. Given these two fields, it is clear that it can be 
difficult at times to associate one particular galaxy with an 
absorption system. One alternative interpretation is that the 
material responsible for the absorption is tidal debris originat- 
ing from both galaxies due to past mergers and harassments 
(similar to the Magellanic stream). Given that the galaxies 
are in close proximity (projected), it is also possible that these 
galaxies share a common gas structure that gives rise to the 
absorption. 

4. 10. Summary II: Disk Halo Model 

In an effort to reproduce the Mgll absorption velocities, we 
used a disk halo model to compute the expected absorption 
velocities. In only one case, Q0229+ 131 Gl, we were able to 
reproduce the full spread of the Mg II absorption velocities in 
a disk halo model. In four other cases, Q0454-220, Q0450- 
132 (G2), and Ql 127-145 (Gl and G2), the velocity region 
with the strongest absorption is consistent with the model. 
However, the halo model of the galaxy G2 toward Q0450- 
132 does not account for roughly 35% of the absorption. In 
the cases of both Q0229+ 131 and Q1127- 145, the models 
cannot account for the unsaturated small cloud structures at 
higher velocities relative to systemic. For each case, the mod- 
els do reproduce some of the absorption velocities, however, 
the disk rotating halo model is insufficient to account for the 
full observed Mg II absorption velocity range. 

We emphasize that the disk halo model applied to the data 
is an extreme version of the spatial-kinematic relationship in 
that (1) all the gas is assumed to rotate at the maximum ob- 
served velocity of the galaxy, and (2) the scale height of the 
models (/z,, = 1000 kpc) is highly unrealistic. These extreme 
conditions were required in order to obtain the greatest degree 
of agreement between the model and the data. If these condi- 
tions are relaxed, the level of agreement we found is dimin- 
ished substantially. None the less, even under these extreme 
and unrealistic model parameters, the disk halo model pro- 
vides insight into the degree at which rotation kinematics can 
account for limited regions of the absorption velocity spread. 

What we learn from the disk halo model is that it is 
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reasonable to suggest that additional dynamical processes 
(such as infall, outflow, supernovae winds, mergers, etc.) 
and/or additional satellites or unidentified galaxies giving rise 
to some of the Mgll absorption contribute to the observed 
velocity spreads. The possibility of unidentified galaxies is 
difficult to quantify, for the assigning of a given galaxy to 
a given absorption system is by its very nature not 100% 
certain. It would, however, require the galaxies we have 
assigned to not be a significant contributor to the absorption 
and/or the unidentified galaxy to have very extended ab- 
sorbing gas. The former statement is not strongly supported 
by a body of previous st u dies (e.g ., Beraeron & Boissl 
1991'; LeBrunetal. 1993; Steidel, Dickinson, & Persso3 
1994; Guillemin & Bergeron 1997; Steidel et al. 1997t 
Churchill, Kacprzak, & Steidel 2005; Tripp & Bowen 2005; 
Zibetti et al.. ,2007; .Chen & Tinker, ,200& .Tinker & Chen 
20081) ' The later statement is based upon the detected galaxies 
in the quasar fields. The statement would not apply to 
putative galaxies below our detection limit that are in close 
projection to the quasar. 

5. ACDM COSMOLOGICAL GALAXY SIMULATIONS 

To further understand the halo gas dynamics and the mech- 
anisms driving the observed gas velocity spread, we investi- 
gate high resolution cosmological simulation of galaxy for- 
mation, which include the dynamical processes of infall, out- 
flow, supernovae winds, mergers, etc. Simulations provide 
the only theoretical means to fully incorporate these dynami- 
cal processes in a cosmological setting. We use the method of 
quasar absorption lines through the simulations to "observe" 
the Mgll absorption kinematics. Here we analyze a single 
z = 0.9 simulated galaxy in detail to study the possible struc- 
tures and mechanisms that give rise to the observed Mg II halo 
gas kinematics. By comparing halo gas kinematics in the sim- 
ulations to the spatial and dynamic processes of the simulated 
galactic environments, we can gain further insights into the 
observed Mgll absorption properties. 

5.1. Description of The Simulations 

The ACDM cosmological simulations are performed us- 
ing the Eulerian Gasdy namics plus N-body A daptive Refine- 
ment Tree (ART) code dKravtsovll 1 9"99l l2003h . Physical pro- 
cesses implemented in the code include star formation and 
stellar feedback, metal enrichment from type 11 and la super- 
novae, self-consistent advection of metals, and metallicity- 
dependent cooling and photoionizat ion heating due to a cos - 
mological ultraviolet background dHaardt & Madaul Il996l) . 
The code does not include radiative transfer, magnetic fields 
or Kelvin-Helmholtz instabilities. The cosmological model 
has a„ = 0.3, 17a = 0.7, h = 0.7, and = 0.9. 

These simulations have high star formation efficiency at the 
resolution scale. This assumption is motivated by studies of 
star formation in simulations of isolated disks with a simi- 
lar resolution (Tasker & Brvan 2006). The gas consumption 
time-scale is ~ 10^ yr for the star-forming cells, however, only 
a small fraction (^1%) of the gas in the disk is forming stars 
in a typical time step of ^40 Myr at z^l. As a result, the 
disk-averaged star formation efficiency is low: the gas con- 
sumption time-sc ale is 10^ yr, consistent with observations 
(|Kennicutt"1998). 

The computational region is a 10 /;"' Mpc co-mo ving box. 
We apply a zooming technique (iKlvpin et al.ll200lT) to select 
a Lagrangian volume of three viral radii centered in a MW- 
size halo at redshift z = 0. The volume is then re-simulated 



from z = 50 with higher resolution and hydrodynamics. The 
high-resolution region has a radius of about 1.5 /;"' co-moving 
Mpc, and has about 4 x 10^ dark matter particles with 5 x 
IO^Mq per particle. The volume is resolved with about 1 .3 x 
10^ hydrodynamic cells with different levels of resolution. 

The combined effect of stellar winds and super- 
nova explosions at the resolved scale dCeverino & KlvpinI 
|2009|) prevent the over-cooling problem of galaxy forma- 
tion at high redshifts (White & Frenk 199 1) and reduce 
the a ngular momentum problem found in early simula- 
tions dNavarro & Steinmetj l2000h . This results in mod- 
els with flat rotation cur ves consistent with observations 
dCeverino & KlypinI l2009l) . and is achieved without typi- 
cal ad-hoc assumptions about the physics at sub-resolution 
scales. These high resolution simulations allow us to resolve 
the regime in which stellar feedback overcomes the radiative 
cooling. By resolving this regime, simulations naturally pro- 
duce galactic scale outflows in star-forming galaxies (Cev- 
erino & Klypin 2010, in preparation) and galaxy formation 
proceeds in a more realistic, although violent way, through 
a combination of cold flows accretion, mergers, and galaxy 
outflows. 

5.2. Simulated Spectra 

To study the Mgll absorption arising in the gas halos within 
the simulations, we employ the following methods. For a 
given gas cell probed by a line of sight through the simula- 
tion box, the total hydrogen density, temperature, and metal- 
licity is used to obtain the Mg II ionization fraction assuming 
photoio nization condi tions. Po st simulation, we use C loudy 
(V96b4. lFeriandll200i ) with the'Haardt & Madau (1996) UV 
background spectrum at the appropriate redshift. The line of 
sight the redshift, Mg II column density, and Doppler b param- 
eter (assuming thermal broadening) are computed for each 
cell. 

Absorption spectra with the instrumental and noise charac- 
teristics of the HIRES spectrograph are generated assuming 
each cell gives rise to a Voigt profile at its line of sight redshift. 
We give each spectrum a signal-to-noise ratio of 50 per pixel, 
which corresponds to a limiting equivalent width detection of 
0.005 A for unresolved lines. The spectrum for each sightline 
is then objectively analyzed for detectable absorption above 
the equivalent width threshold of 0.02 A, which corresponds 
to log A^(Mg II) = 1 1 .7 cm"^ for b = 5 km s"' . The mean opti- 
cal depths (mean redshifts), rest-frame equivalent widths and 
velocity widths, and other quantities are then measured (see 
IChurchill & Vogti,2001) . The velocity zero point of the simu- 
lated absorption lines is set to the line of sight velocity of the 
simulated galaxy (center of mass of the stars). 

To examine the 3D spatial and kinematic properties of gas 
giving rise to Mg II absorption, we identify Mg II "absorbing 
gas cells" along each sightline as those which contribute to 
detected absorption in the simulated spectra; they are defined 
as cells that align within the range of line of sight velocities 
of the a bsorption. We account f or multiple kinematics sub- 
systems (Churchill & \^l20()Tl) . regions of absorption sepa- 
rated by continuum. 

5.3. Discussion of Simulated Galaxy Observations 

We focus on a single typical galaxy at z = 0.923. The galaxy 
star formation rate is SFR = 3.5 M0 yr"'. The galaxy has 
a maximum rotation velocity of 180 km s"' when observed 
edge on. Based upon the the Tully-Fisher relation, we derive 
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Fig. 9. — The integrated total hydrogen column density, N(H), over a 220 kpc cube is shown for an z = 0.923 simulated galaxy viewed edge on. The direction 
of the simulated quasar lines of sight are perpendicular to the plane of the image. Squares of increasing size are plotted where Mgll absorption was detected 
along the line of sight in the simulated quasar spectra. We apply an equivalent width detection limit of Wr(2796) > 0.02 A. The four square sizes indicate, in 
increasing order, Mgll absorption equivalent width bins of; 0.02 < W,(2796) < 0.3, 0.3 < lVr(2796) < 0.6, 0.6 < Wr(2796) < 1.0, and 1.0 < W,(2796) < 3.0, 
respectively. 




Fig. 10. — Same as Figure|9] except the integrated neutral Hi column density, N(Hl), is shown. The direction of the simulated quasar lines of sight are 
perpendicular to the plane of the image. Note that small DLA regions, having N(Hl) > 10-"^ cm"^, are seen beyond ~ 85 kpc. 
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Fig. 1 1 . — Same as Figure|9]for the same galaxy, except viewed face on. The direction of the simulated quasar lines of sight are perpendicular to the plane of 
the image. 




Fig. 12. — Same as Figure flOl for the same galaxy, except viewed face on. The direction of the simulated quasar lines of sight are perpendicular to the plane of 
the image. 



18 



Kacprzak et al. 






Tr,-,Kr,-|v,irr,-,-r 


,T7-,-;i7,-r|-,-;7rr,-,-| 


TT, -. I , , 

















Fig. 13. — A 45 X 45 kpc region extracted from the lower left quadrant of Figure 1 101 The region shows two satellite galaxies and their tidal streams. The 
velocity zero point of the absorption profiles is the galaxy systemic velocity. We enforced a detection sensitivity limit of W,-(2796) > 0.02 A. The majority of the 
Mgll absorption arising in this tidal stream is radially infalling towards the galaxy. 



a luminosity of Lg = OALg. The average luminosity for our 
sample is {Lb) = 0.58Lg excluding the 8Lg galaxy. Thus, this 
simulated galaxy is well representative of our observational 
sample. 

The simulated galaxy is probed with with a square grid of 
sightlines at intervals of 7.5 kpc that span -110 kpc to +1 10 
kpc on the "sky" and 14 Mpc along the line of sight^. There 
are 900 total sightlines. The highest resolution of the adaptive 
mesh at z = 0.923 is 225 pc. The gas contributing to detectable 
Mg II absorption is found in a range of cell resolutions from 
225- 1815 pc with the majority of the gas arising in cells of 
resolution 905 pc. We examine simulated quasar lines of sight 
for this galaxy with three different inclinations, / = 5° (face- 
on), 45°, and 85° (edge-on). 

In Figure |9l we show the integrated total hydrogen column 
density, N(H), over a 220 kpc cube for the edge-on view of 
the galaxy. There is no absorption detected along the line of 
sight outside of this cube. The galaxy is clearly not in iso- 
lation. In the image, filaments and tidal stream material can 
be seen. Several low mass satellites galaxies are also in the 
process of interacting with the main galaxy. We have super- 
imposed squares over the sightlines where Mgll absorption 
was detected in simulated quasar spectra. Increasing square 
sizes indicate the absorption strength in four bins: 0.02 < 
W,(2796) < 0.3, 0.3 < W,(2796) < 0.6, 0.6 < W,(2796) < 1 .0, 
and 1.0 < Wr(2796) < 3.0. Out of the 900 lines of sight, Mgll 
absorption was detected in 87 for the edge-on case. Note 

* For this and all subsequent discussions, spatial quantities are quoted as 
proper lengths 



that the absorption is not distributed ubiquitously on the sky 
around the galaxy, but traces the various structures around the 
galaxy. The covering fraction of the Mg II gas is low 10%). 

In Figure [TOl we show the integrated neutral hydrogen col- 
umn density, N(Hl), over the same cube. The N(Hl) range 
from roughly 10'^ to 10^^ cm"^. The structures that are as- 
sociated with absorption stand out more in Hi. We see as 
expected dChurchill et alj|1999 ) that Mgll absorption is de- 
tected in regions with N(Hl)> lO"''^ cm"^. Regions in which 
N(Hl)> 10"'-^ cm"^ that fall between the line of sight grid 
sampling are also expected to produced Mgll absorption. 
Note that as one goes to lower N(Hl), the gas covering frac- 
tion increases. 

In Figures [TTlfT2l we present N(H) and N(Hl), respectively, 
for the face-on view of the galaxy. Out of the 900 lines of 
sight, Mgll absorption was detected in 96. The covering frac- 
tion remains roughly 10%. Again, the absorption primarily 
arises in streams and filaments. 

We do not show the 45 degree view of the galaxy here. Out 
of the 900 lines of sight, Mg II absorption was detected in 124, 
and the covering fraction is roughly 20%. 

Overall, there are many complicated structures that reside 
within what we classically call a halo. Note that we find pock- 
ets of DLA Hi column densities [N(Hl)> cm"^] out to 
^100 kpc. Although, the covering fractions of these dense re- 
gions are low, their presence suggests that it is possible to ob- 
serve DLA absorption at high impact parameters. These large 
impact parameter DLA systems arise from the inner regions 
of satellite galaxies. These satellites have a luminosity range 
0.003L* < Lb < 0.006L*. Assuming standard K-corrections, 
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Fig. 14. — A 37.5 x 45 kpc inner central region of the edge-on galaxy seen in Figure fTol The velocity zero point of the absorption profiles is the galaxy systemic 
velocity. We enforced a detection sensitivity limit of W,-(2796) > 0.02 A. In the inner ~ 15 kpc of the galaxy center we find some gas outflowing at ~ 200 km s"' . 
The Mgll absorption profiles produced by these outflows are similar to those seen for Q0002+051 Gl (Figure|2})) and Q0836+ 113 Gl (Figure|4l)). 



their apparent magnitude in the HST F702W filter range be- 
tween 28 < infji)2w ^ 29. These DLA producing satellites 
are well below the typical detectability of WFPC-2/HST and 
ACS/HST. Also note that the satellite in the upper right cor- 
ner of Figure [12] has an Hi morphology similar to galaxies 
seen on their first pass through a cluster, where gas is being 
stripped (e.g.J Chung et al.ll2007h . 

Focusing on the two small satellite galaxies and tidal stream 
in a 45 X 45 kpc region of the lower left quadrant of Fig- 
ure [TOl we present an example of simulated Mgll A2796 ab- 
sorption spectra in Figure [13] The velocity zero point of 
the simulated spectra is set to the galaxy systemic velocity. 
The absorption profiles are quite similar to those detected in 
Mgll surveys. The simulated profiles have velocity spreads of 
~ 25-100 km s~', some comprising multiple kinematic com- 
ponents and complex structures. The Mgll absorption spans 
across the galaxy systemic velocity or resides blueward of the 
galaxy systemic. Absorption arising near the galaxy systemic 
velocity or entirely to one side is consistent with the majority 
of our observational data. 

In Figure [14] we show the Mgll A2796 absorption spectra 
over the 37.5 x 45 kpc central region of the edge-on galaxy 
shown in Figure [TO] Individual absorption profile velocities, 
in the inner regions of the galaxy, show a variety of complex 
kinematics and optical depths. In our observational sample 
of galaxies, we have no galaxy-absorber systems that have 
impact parameters less than 25 kpc. In the simulations, the 
inner ~ 15 kpc contains some ^ 200 km s"' outflowing gas. 
The outflows are not strong, but their signatures are reflected 
by the complex kinematics of the absorption profiles in the 



inner regions. The reason why there is a line of sight near the 
center of the galaxy that does not produce Mg II absorption is 
because it is dominated by 10^ K gas. The profiles reduce in 
kinematic complexity rapidly as impact parameter increases. 

To the degree that the simulations reflect the reality of the 
gaseous environments around galaxies, we find that Mg II ab- 
sorption arises in many types of structures (metal enriched fil- 
aments, minor satellites, tidal streams, and within the region 
of the galaxy itself). As inferred from our simulated absorp- 
tion line survey, and guided by analysis of the 3D simula- 
tions (examples provided by Figures l9i-[T4]i. the simulations 
are not suggestive of Mg II absorption arising from spherical 
or disk-like halos with uniform gas covering fraction. We re- 
visit the spatial and kinematic distribution of the absorbing 
gas in § 15.51 

5.4. Disk Halo Models of The Simulated Galaxy 

Given the structures shown in Figures l9i-[T2] it would seem 
to be unrealistic to treat the halo as a monolithic thick disk. 
However, we have modeled all the simulated Mg II absorption 
profiles for each of the three galaxy inclinations with the disk 
halo model (§|4]l. 

We find that for the edge-on view, out of 87 lines of sight 
with detected Mgll absorption, 45% have kinematics consis- 
tent with the model (the full range of velocities can be ex- 
plained). For / = 45°, out of 124 lines of sight, 26% have 
kinematics consistent with the model. For the face-on view, 
out of 96 lines of sight, only one is consistent with the model. 
This is not surprising, since the projected maximum rotation 
velocity is small compared to the absorption spreads. 
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Fig. 15. — The spatial distribution of N(MgII)< 1 1.5 cm"^ gas contributing to the Mgll absorption along the lines of sight. The edge-on galaxy is located 
at the origin; the black horizontal line represents a line of sight passing through the galaxy center. The observer is looking along the positive x direction. All 
simulated lines of sight are parallel to the .v-axis. The absorbing gas is color coded as a function of line of sight velocity relative to the galaxy systemic velocity, 
as coded in the legend. Red absorbing gas is moving away from the observer; blue is moving toward. 




Fig. 16. — Same as Figure [15] except the velocity color coding is for the radial velocity component with respect to the galaxy. Red absorbing gas is outflow; 
blue is inflow. The majority of the Mg II absorption arises in filaments and tidal streams and is infalling towards the galaxy. 
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Fig. 17. — Velocity distributions of Mgll abso rbin g g as c ells within the 
220 kpc cube of the simulation shown in Figures [Ts] and [T6l Gas within a 
20 kpc radius of the galaxy is omitted. — (a) The total velocity distribution. 
— (b) The radial velocity distiibution. — (c) The distribution of the ratio of 
radial to total velocity. A ratio of v,- / Viot = it 1 indicates pure radial motion. 



The conclusion to be drawn here is, that even if the halo gas 
does not rotate as a monolithic disk in the simulations, Mg ll 
absorption detected along some lines of sight can still appear 
to be consistent with a disk halo model. 

5.5. Halo Gas Spatial and Velocity Distribution 

To understand what kinematic mechanisms are responsible 
for the Mg II absorption velocity spreads measured in the sim- 
ulated spectra, we examine the velocities and spatial distribu- 
tions Mgll absorbing gas. 

In Figure [15] we present the 3D spatial distribution of 
N(Mgll)> 10"-^ cm"^ gas contributing to the Mgll absorp- 
tion along the lines of sight. The edge-on galaxy is located 
at the origin and the lines of sight are parallel to the x- 
axis. Gas with N(Mgll)< 10"^ cm"^ is not shown for clarity 
(which corresponds to an equivalent width limit of 0.012 A for 



b = 5 km s"'). From this view, the tidal streams and filaments 
can be visually discerned. There are also isolated clouds that 
produce Mg II absorption. Regions that contribute to absorp- 
tion have physical sizes of 225 pc (simulation resolution) to 
~ 20 kpc. All the absorption along the lines of sight occurs 
within ±110 kpc from the center of the galaxy. 

The gas is color coded as a function of line of sight velocity 
relative to the galaxy systemic velocity. Gas that is colored 
red (blue) is moving away from (towards) the observer. The 
gas has line of sight velocities ranging from -160 < v/o, < 
+160 km s"' . The line of sight velocity dispersion in the inner 
regions near the galaxy appears inconsistent with disk rota- 
tion. The somewhat randomized velocities reflect the winds. 
Upon carefully examining the gas velocities along a particular 
line of sight, velocity gradients can be observed. For exam- 
ple, the region seen in Figure [T5]at x = 100, y = -50 kpc, has 
velocity gradients of ~ 30-40 km s"' . 

In Figure [16] we present the same three dimension spatial 
distribution shown in Figure [15] except that the gas is color 
coded as a function of radial velocity relative to the galaxy. 
Gas colored red (blue) has a radial velocity component that 
is outflowing from (infalling towards) the galaxy. This com- 
bined spatial and kinematic representation provides an holis- 
tic view of the halo dynamics. The gas has radial velocities 
ranging from -200 < < +100 km s"' . If one focuses on the 
filament structure (x = kpc, y = - 100 kpc) in the plane of 
the galaxy, one can see the strong velocity gradient. The gas 
along the filament far from the galaxy increases in infall ve- 
locity from 70 to ~ 200 km s"' as it approaches the galaxy 
center The same dynamics can be seen along the tidal stream 
(jc = 0-50 kpc, y = -50 kpc) originating from the two small 
satellite galaxies. 

In Figure WTh . we show the probability distribution of the 
total gas velocity of the Mgll absorbing gas cells. What we 
call the probability distribution is the area normalized fre- 
quency distribution that we detected in the simulations. The 
total velocity is the magnitude of the gas velocity vector rel- 
ative to the galaxy. Since our observational data have only 
impact parameters greater than ^ 20 kpc, we exclude all ab- 
sorbing gas within D < 20 kpc of the simulated galaxy. The 
velocities range from 50 < v,o/ < 250 km s~' with a peak at 
185 km s"' and secondary maxima at 100 kms"'. In Fig- 
ure [TtJ?, we show the probability distribution of radial veloc- 
ities of the Mgll absorbing cells. The velocities range from 
-200 < vv < +1 10 km s"' with a maximum at -160 km s"'. It 
appears that, beyond D > 20 kpc, most of the gas is infalling 
towards the galaxy and very little is outflowing. In Figure[T7b. 
we show the ratio of the radial to the total velocity. The bulk 
of the gas is dominated by radial infalling velocities. 

Drawing from the number of lines of sight we have through 
the simulations, we produced the probability distribution of 
absorption velocity offsets from the galaxy systemic veloc- 
ity (Av = Vahs-Vgai kms"') using the simulated absorption 
profiles. The quantity Vats is obtained by calculating the 
optical depth weighted mean of the profiles (the velocity at 
which there is equal optical depth to both sides along the 
profiles). In Figures [TSl;./?. and c, we show the velocity 
offset probability distributions for the edge-on, / = 45°, and 
face-on orientations, respectively. For the edge-on case, the 
absorption velocity offset ranges from ±100 km s~' with a 
strong peak around +30 km s"'; it is highly probable to de- 
tect absorption to one side of the galaxy systemic velocity. 
For the / = 45° case, the velocity spread increases to about 
±200 kms"' and develops multiple peaks at ±150 kms"' 
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Fig. 18. — Probability distributions of the mean velocities of the simu- 
lated Mgll absorption profiles arising in the lines of sight, for three galaxy 
orientations. The quantity v„i,s is the optical depth weighted mean of the pro- 
files. Lines of sight with impact parameters less that 20 kpc are omitted. — 
(a) (' = 85°, edge-on. — (b) i = 45°. — (c) = 5°, face-on. The absorption 
profiles span the rotation velocities of the simulated galaxy (iv = 180 km s"' ). 

and at ^ +50 km s~' . The probability of detecting gas at the 
galaxy systemic velocity is significantly smaller than over the 
range |50- 150| km s"' . For the face-on case, the velocity dis- 
persion is still around ±200 km s"' . Again, it is unlikely that 
absorption will be detected at the galaxy systemic velocity. 

Since most of the Mgll gas arises between ±200 km s~', 
with very little at the galaxy systemic velocity, and since 
most projected galaxy rotation curves have maximum ve- 
locities of ±200 kms"', it may not be a surprise to ob- 
serve Mgll absorption aligned with the observed galaxy ro- 
tat ion curve. These re sults are consistent with the findings 
of iBouche et al.l (l2007h . who detected galaxy Ha emission 
within ±200 km s"' of the optical depth weighted mean Mg 11 
absorption redshift. 

5.6. Shortcoming of the Simulations 

The technique of quasar absorption lines through cosmo- 
logical simulations is one of several promising approaches to 
understanding the dynamics of galaxy halos. At the present 
time, simulations of galaxy formation in the cosmological 
context still need to achieve greater accuracy for modeling 
stellar feedback. For this study, we employed a feedback 



recipe that successfully results in extended metal enriched 
gas around galaxies. These simulations result in an equiv- 
alent width distribution with an under abundance of larger 
equivalent widths and a relative over abundance of smaller 
equivalent widths. They also under predict the observed 
Mgll mean absorption covering fra ction range of 0.2-0.9 
(Tripp & Bowen 2005t [Chen & Tin ker 2008; Kac przak et all 
2008; Barton & Cookell2009h . The covering fractions for all 
absorption above 0.02 A are as follows: (1) / = 90° (edge- 
on): total 10%, weak 6%, strong 3%; (2) / = 0° (face-on): 
total 10%, weak 6%, strong 3% (3) / = 45°: total 14%, weak 
10%, strong 5%. 

These mismatches with observations could either be a result 
of the method in which Mg 11 column densities are determined 
in the simulations or observational biases. In the simulations, 
the determination of the Mg II ionization fraction may be un- 
derestimated due to the fact that we do not account for shield- 
ing of UV photons in the ionization corrections. It is also 
possible that the resolution of the simulations may influence 
the derived Mgll column densities and that higher resolution 
may in fact lead to higher column densities. We aim to an- 
alyze such issues in future work. A possible observational 
biases that galaxies are selected, identified, and assigned to 
already known absorption sys tems. This may elevat e the in- 
ferred covering fraction (e.g.. lTripp & Bowenll2005h . We do 
emphasize, however, that in the simulations we do detect 
Mgll absor ption out to ~ 120 kpc, as seen in current ob- 
servations (Churchill. Kac przak. & SteideJ2005tlZibetti et al.l 
2007; Kacprzak et al. 20081). 

As an additional caveat, we also remark that the experiment 
to examine the spatial and kinematic relationship between the 
galaxy and the Mgll absorbing structures in the simulations 
is very different than the observational experiment in one re- 
gard. We examine a grid of sight lines through a single simu- 
lated galaxy, which is in a unique environment and undergoes 
a unique evolution in the IGM. The observational data, on 
the other hand, sample a single line of sight through various 
galaxies in various environments and with various evolution- 
ary histories and with random orientation through the galaxy 
and environment. 

We re-emphasize that the simulations, as we have applied 
them here, provide a first view of the types if physical struc- 
tures and their spatial and kinematic relationship to a galaxy in 
the cosmological environment that give rise to Mg II absorp- 
tion. Though the details of the Mg II absorbing properties are 
not yet tuned to observations (and we mention that our experi- 
ment is only for a single galaxy using an incomplete sampling 
of its halo- a full comparison with survey statistics is not en- 
tirely applicable), the dynamical structures themselves in the 
simulations are robust. That we find some parts of these struc- 
tures give rise to Mg II absorption is not outside expectation 
given the densities and temperatures of the gas. The upshot is 
that, given these discrepancies and concerns, the simulations 
are not expected to provide a fully quantitative comparison 
with the data. However, the simulations do provide a fully 
self-consistent galaxy model from which we acquired valu- 
able insights for interpreting Mgll absorption line observa- 
tions. 

5.7. Summary III: Galaxy Simulations and Halo Gas 
Distributions 

In our study of z = 0.923 galaxy at three inclination an- 
gles with respect to the simulated quasar lines of sight, we 
have detected Mgll absorption in a variety of structures. 
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Mgll absorption was detected in inflowing metal enriched 
filaments, tidal streams, small satellites, and gas around the 
host galaxy. The absorption resides in a "halo" of about 
100 kpc in size. The types of structures that form in the 
simulations (see Figures l9UT2li are a challenge to models 
in which t he absorption is assigned to thick disks (e.g., 
ICharlton"& Churchill 19 9^ and symm etrically distributed 
halos (e.g.. Mo & Mirald a-Escudel ll996). The structures also 
may provide further guidance for expanding upon radial den- 
sity d ependent halo occupation models (e.g. iTinker & ChenI 
l2008h . The spatial distribution of filaments, tidal streams, and 
satellite galaxies are asymmetric, patchy (low volume filling 
factor), complex, and part of a cosmological setting. 

DLA Hi column densities are seen out to ^ 100 kpc, and 
arise in low mass satellite galaxies. These galaxies are below 
the detection limits of even deep HST images. Although, the 
covering fractions of these dense regions are low, this might 
explain why some bright galaxies at DLA redshifts are found 
at large impa ct parameters (e.g., the 3C 336 z = 0.656 DLA at 
D = 100 kpc: ISteidel et"anil997h . 

In the simulations, the kinematics are closely coupled to the 
gas structures. Metal enriched tidal streams and filaments are 
dominated by infall to the central galaxy, and these structures 
are selected by Mgll absorption (see Figures [TS] and [T6]l. In 
fact, we find that gas giving rise to the Mgll absorption is 
dominated by inflow (see Figure [TT] ). Inflow velocity gra- 
dients are apparent, such that the infall increases as the gas 
approaches the galaxy. This inflow is not spatially symmetric. 

Similar to the observations, the velocity and spatial distri- 
butions from the simulations conspire to give rise to Mgll 
absorption to one side of the galaxy systemic velocity even 
though the absorbing gas is not rotating with the star forming 
component (see Figures[T3]andfT4l). The velocity offset prob- 
ability distribution (relative to the simulated galaxy) spans 
^ 200 km s"' with lowest probability of detecting Mg II at 
the galaxy systematic velocity (see Figure [TS]). Thus, the fact 
that we observe Mgll absorption velocities consistent with the 
galaxy rotation curves may be a natural consequence of the 
spatial and kinematic distributions of gas in complex environ- 
ments surrounding galaxies. 

The current feedback recipes used in the simulations suc- 
cessfully produces in extended metal enriched halo gas 
around galaxies. However, for this single galaxy we de- 
rive covering fractions that are low er than the current ob- 
served mean values of 0.2 - 0.8 ( Tripp & Bowen 2005; 
Chen & Tinkedl2008t iKacprzak et al. 2008; Barton & CookQ 
2009^ . In this particular galaxy halo, we also find an un- 
der abundance of larger equivalent widths and a relative over 
abundance of smaller equivalent widths. We have yet to de- 
termine if some of these mismatches are due to observational 
biases in current studies or issues with the feedback/baryon 
physics implemented in the simulations. However, the struc- 
tures in the simulations still provide valuable insight for inter- 
preting Mgll absorption line observations. 

6. LARGE scale GALACTIC WINDS: OBSERVATIONS AND 
SIMULATIONS 

In our analysis of the simulated galaxy, we do not find 
systematic rotation near the galaxy, as was suggested by 
iLanzetta & Boweiil (|1992|) . but that the gas kinematics reflects 
processes such as winds and chimneys. These are characteris- 
tic of stellar feedback processed that successfully circumvent 
the over cooling and ang ular momentuni prob lems that have 
plagued simulations CCeverino & Klvpinll2009l) . 



Observationallv.' Weiner et aTl (l2009l) has inferred outflows 
in Mg II absorption associated with \ star-forming galax- 
ies (10 < SFR < 40 M0 yr"'). The outflow velocities are 
proportional to SFR^^ ^. The star formation rate in our sim- 
ulated galaxy is 3.5 M0 yr"'. Calibrating to their obser- 
vations, our galaxy is expected to have winds in the range 
220- 350 km s"' . Our simulated Mgll profiles indicated some 
outflowing gas with velocities of ~ 200 km s"' in the inner 
15 kpc (independent of galaxy disk inclination). The absorp- 
tion profiles produced by the outflows are saturated and span 
both sides of the systemic velocity of the simulated galaxy 
(see Figure O. Thus, the saturated Mgll absorption profile 
that spans both sides of the galaxy systemic velocity might be 
direct signature of outflows. 

Two galaxies in our sample, Q0002 + 051 Gl and Q0836 + 
113 Gl, have saturated profiles that span both sides of the 
galaxy systemic velocity. These saturated, symmetric pro- 
files are also observed at high redshift (z ^ 2) and span the 
galaxy systemic velocity. The absorption is attributed to large 
scale galactic outflows which can be detected at least out to 
impact parameters of ^ 70 kpc (Steidel et al. 2010, in prepa- 
ration). The outflows are determined to be more or less sym- 
metric and radial, giving rise to symmetric absorption pro- 
files. These high redshift galaxies do not exhibit substantial 
kinematic structure and are instead consisten t with dispersion - 
dominated velocity fields around 80 km s"' (iLaw et al.ll2007h . 
similar to the shear observed for Q0002+05 1 Gl and Q0836+ 
1 1 3 G 1 . Given the similarities between these two galaxies and 
those at high redshift observed to have outflows, one possibil- 
ity is that large scale outflows are responsible for the absorp- 
tion velocities associated with both galaxies. To explore the 
outflow scenario for Q0002+051 Gl and Q0836+ 1 13 Gl, we 
examine the galaxy star formation rates. 

For Q0836+ 113 Gl, we measure a [Oil] line flux 
of 3.4 ± .43 X lO""* e rg s"| cm'^. We use the rela- 
tion from iKewley et al.l (l2004t) to obtain a star formation 
rate of 6.5 M^yrK For Q0002 + 051 Gl we determined 
the star formation from the UV flux at 1700 A which 
was d erived from the j ?-band magnitude of the galaxy 
(Steid el. Dickinson. & Per sson 199 4j). The r obustness of this 
method has been demonstrated by .Erb et alJ (120061 . From the 
measured UV flux of 2.3 ±0.5 x 10"^^ erg s~' cm"^, we de- 
rive a SFR of 6.5 MQyr"' . The SFRs are not corrected for the 
internal extinction of the galaxies and are thus lower limits. 

Using the results of IWeiner et alJ (120091) . both galaxies 
are expected to have outflows with velocities in the rage 
270-430 kms"'. Both galaxies are at impact parameters 
of ~ 25 kpc, which is near the boundary of where we lose 
the signature of outflows for our simulated galaxy. Although 
Q0836+ 113 Gl and Q0002 + 051 Gl have Mgll absorption 
velocity widths of 290 and 470 km s"' , respectively, which 
are compara ble to the outflow ve locity range predicted from 
the results of'Weiner et al.' (l2009l) . 

Jieckman (2002, 2003) discusses that the star formation 
per unit area is a more relevant indicator of galaxy out- 
flows. These outflows are ubiquitous in galaxies where 
the global star-formation rate per unit area exceeds E = 
0.1 Mq yr"' kpc"^, where the area is defined by the half 
light radius of the galaxy. This criteria applies to local star- 
bursts and even high redshift Lyman Break galaxies. The 
ISM entrained in the winds have outflow speeds of ^ 100 
to - 1000 km s"'. For Q0002 + 051 Gl we obtain a E > 
0.35 Mq yr-i kpc'l For Q0836+ 113 Gl we obtain a S > 
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0.14 M0 yr~' kpc~^. Thus, both of the galaxies are expected 
to have outflow signatures. 

We estimate that the simulated galaxy has a E = 
0.08 Mq yr"' kpc"^, which is slightly less that the criterion 
of Heckman. This is consistent with our outflow velocities 
derived from the Mgll absorption profiles since we do not 
see evidence of strong large scale outflows. This particu- 
lar simulated galaxy may not be well representative of the 
Q0002 + 051G1 andQ0836+113Gl galaxies. 

In summary, galaxies Q0002 + 051 Gl and Q0836+ 113 Gl 
are kinematically similar to high redshift absorption selected 
galaxies. The SFRs and Es for both galaxies exceed the limits 
were strong outflows are expected. Given the large impact 
parameters that outflows are detected at high redshift, it is 
quite possible the observed Mg II absorption kinematics for 
galaxies Q0002 + 051 Gl and Q0836+ 113 Gl are signatures 
of outflowing gas. 

7. CONCLUSIONS 

We have examined and compared the detailed galaxy and 
Mgll absorbing kinematics for a sample of 10 intermedi- 
ate redshift galaxies. The galaxies have a wide range of in- 
clinations and orientations with respect to the background 
quasar. The galaxy-quasar impact parameters range from 
26 <D < 108 kpc. The galaxy rotation curves were obtained 
from ESI/Keck spectra and the Mg II absorption profiles were 
obtained from HIRES/Keck and UVES/VLT quasar spectra. 
In an effort to compare the relative kinematics, we used a thick 
disk halo model to compute the expected absorption velocities 
through a monolithic gaseous halo. 

To obtain theoretical insights into the gas dynamics and spa- 
tial distribution of halos, we used the technique of quasar ab- 
sorption lines to analyze Mg II absorption around a galaxy in a 
high resolution cosmological simulation of galaxy formation. 
The galaxy was probed with a square grid of sightlines at in- 
tervals of 7.5 kpc that span -1 10 kpc to +1 10 kpc for a total 
of 900 sightlines. We examined this galaxy at three different 
inclinations, face-on, 45°, and edge-on. 

Our mains results can be summarized as follows: 

1 . For all ten galaxies, the velocity of the strongest Mg II 
absorption component lies in the range of the observed 
galaxy rotation curve. In seven of ten cases, the Mg II 
and Mg I absorption velocities reside fully to one side 
of the galaxy systemic velocity. The strongest absorp- 
tion usually aligns with one arm of the rotation curve. 
In the three remaining cases, the absorption velocities 
span both sides of the galaxy systemic velocity. Two 
of those three (Q0002 + 051 Gl and Q0836+ 113 Gl) 
have strong saturated absorption on both sides of the 
galaxy systemic velocity. The third (Q1127- 145 G3), 
has two very weak clouds, and therefore probes low 
column density gas. 

2. For galaxies Q0002 + 051 Gl and Q0836+ 113 Gl, we 
have determined that large scale galactic outflows might 
be giving rise to the observed Mgll absorption kine- 
matics. Both galaxies have a SFR 6.5 Mq yr"' and 
E > 0.35 Mq yr-i kpc'^ and E > 0.14 Mq yr"' kpc'^, 
respectively. These SFRs and Es are typically found for 
galaxies exhibiting outflow velocities of several hun- 
dred km s"' . The Mg II absorption velocities associated 
with the two galaxies span both sides of their systemic 
velocity. Such profiles have been interpreted, in both 



our simulations and at high redshift, as signatures of 
outflows. 

3. We find that the observed Mgll absorption velocity 
spread and optical depth distribution may be a func- 
tion of galaxy inclination. Galaxies with higher in- 
clination exhibit a Mgll absorption velocity spread of 
^ 300 km s"' with a somewhat even distribution of op- 
tical depths, whereas, galaxies with lower inclinations 
exhibit a narrower velocity spread of 100 km s~' with 
a clear optical depth peak at v ^ 100 km s"' . A K-S test 
shows that the Mg II optical depth distributions for the 
high and low inclination bins are not consistent at the 
3.3(7 level. These results suggest that the absorbing gas 
is either disk-like or the spatial distribution and kine- 
matics of the structures producing the absorption (i.e., 
filaments, tidal streams, satellites etc.) are closely cou- 
pled to the disk orientation. 

4. We employed simple rotating disk halo models to ex- 
amine whether disk-like rotation is consistent with the 
observed galaxy-gas kinematics. For model parame- 
ters that allow for a 1 Mpc gas scale height and max- 
imum rotation velocity (rigid rotation) the the bulk of 
the observed absorption kinematics can be explained 
by co-rotation with the galaxy. In all cases, the ro- 
tating disk halo models we present are unable repro- 
duce the full spread of observed Mgll absorption veloc- 
ities. This model is a highly unrealistic representation 
of galaxy gas. When the parameters are relaxed to bet- 
ter reflect reasonable gas scale heights and a slowing of 
the rotation speed with height above the disk plane, the 
relative proportion of the gas velocity spread that can 
be made consistent with galaxy co-rotation diminishes 
such that some absorbers cannot have but a tiny frac- 
tion explained by co-rotation. In this simple scenario, 
even if some of the absorbing gas arises in a thick disk, 
what we learn from the exercise is that some additional 
type of dynamical process (such as infall, outflow, su- 
pernovae winds, etc.) must be invoked to explain the 
range of absorption velocities hat cannot be made con- 
sistent with the simple rotating disk halo model. 

5. In two quasar fields, we find pairs of galaxies that align 
in velocity within 100 km s"' of a single, saturated 
Mgll absorption system. For one case, the observed 
velocity range of the strong saturated component can be 
explained by a rotating disk model only if both galaxies 
contribute to the absorption. This challenges the idea 
that an individual Mg II absorber can be assigned to a 
single galaxy, and understood as an isolated halo. 

6. In the simulations, Mgll absorption selects gas struc- 
tures such as metal enriched tidal streams, filaments, 
small satellite galaxies, and the region within ^ 20 kpc 
of the galaxy. Together, these structures extend roughly 
^100 kpc around the galaxy, suggesting that galaxy 
"halos" are a complex composite of the these various 
structures. 

7. For this simulated galaxy the Mgll covering fraction 
is ^ 10%, which is below th e current observational 
estimated means of 2 0-80% dTj ipp & Bowed 120051: 
iKacprzak et al] l2008t iBarton & Cooke 2009^ This 
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may reflect a need for additional tuning of the feed- 
back/baryon physics in the simulations, or indicates 
current observational biases. 

8. In the simulations, DLA Hi column densities arise in 
low mass satellite galaxies at impact parameters as large 
as ^ 100 kpc. These galaxies are below the detection 
limits of deep HST images. Although, the covering 
fractions of these dense regions are low, this might ex- 
plain why some bright galaxies at DLA redshifts are 
found at large impact parameters. 

9. In the simulations, the majority of the Mgll absorb- 
ing gas is infalling in filaments and tidal streams to- 
wards the galaxy with velocities between -200 < v,. < 
-180 kms"'. The velocity offset probability dis- 
tribution (relative to the simulated galaxy) spans ^ 
200 km s~' with lowest probability of detecting Mgll 
at the galaxy systematic velocity. Thus, observed Mg II 
absorption velocities can fall within the range of the 
galaxy rotation curve velocities, even though the gas 
arises in a variety of kinematics structures. 

The gas structures selected by Mg II in the simulations (see 
Figures l9l-[T2li cannot be described as simple thick disks or 
spherical halos. If the simulations reflect reality, it would ap- 
pear that Mg II absorption arises in large ~ 100 kpc halos that 
are built from the local cosmological environment of a mod- 
erate mass galaxy. Complicating the picture is the fact that we 
find groups and pairs of galaxies that align in velocity within 
~ 100 km s"' of a single Mgll absorption system. This chal- 
lenges the idea that an individual Mgll absorber can be as- 
signed to a single galaxy or understood as an isolated halo. 
Though considered subcomponents of halos, smaller scale 
structures like the Magellanic-type galaxies and tidal streams, 
may c ontribute significantly to the detections of Mg II absorp- 
tion (%^itan[Sl6'; 'Kacprziketin|200^ These consid- 
erations lead us to suggest that galaxies and Mg II absorbers 
should be studied and modeled in a environmental context if 
they are to be fully understood. 

In the simulations, the kinematics are closely coupled to 
the gas structures (i.e., filaments, tidal streams, small satel- 
lite galaxies, and the inner 20 kpc of the central galaxy). As 
observed in our data, the simulated Mg II absorption veloci- 
ties fall within the range of the galaxy rotation velocities, and 
rarely at the galaxy systematic velocity. Thus, the simulations 
suggest that observing Mg II absorption velocities consistent 
with the galaxy rotation curves can naturally occur even if the 



absorption arises in many different structures in the complex 
environment of the galaxy. It is these structures that comprise 
halos. 

A natural extension of the work presented here would 
be to perform a similar study (simulations and observa- 
tions) that incorporates the kinematics of higher ionization 
Civ AA1548,1550 and Ovi AA1031,1037 doublet absorp- 
tion. These ions probe lower density and/or higher tempera- 
ture structures and provide a more comprehensive view of the 
gaseous environment around galaxies. Future observations 
with the Cosmic Origins Spectrograph are perfectly suited for 
the galaxy sample presented in this paper It is also important 
to expand the number of galaxy environments studied in the 
simulations. 
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